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SECTION  1 


INTRODUCTION 


Multiobjective  optimization  has  recently  been  acknowledged  as  an  advanced 
design  technique  in  structural  optimization.  Most  of  the  real  world  problems  are 
multidisciplinary  and  complex,  there  are  always  more  than  one  important  objec¬ 
tive  function  in  the  problem.  In  order  to  accomodate  many  conflicting  design 
goals,  one  needs  to  formulate  the  optimization  problem  with  multiple  objectives. 
Formulation  of  a  structural  optimization  problem  consists  of  constructing  a  math¬ 
ematical  model  that  describes  the  behavior  of  a  physical  system  encompassing 
the  problem  area.  The  main  task  is  the  determination  of  the  choice  of  the  design 
variables,  objectives  and  the  constraints.  Sometimes  only  one  dominating  crite¬ 
rion  may  be  a  sufficient  objective  for  minimization.  This  is  true  especially  if  the 
other  requirements  can  be  presented  by  equality  amd  inequality  constraints.  But 
generally  the  choice  of  the  constraint  limits  may  be  a  difficult  task  in  a  practical 
design  problem.  These  allowable  values  can  be  rather  fuzzy  even  for  common 
quantities  such  as  displacements,  stresses,  and  natural  frequencies.  If  the  limit 
cannot  be  determined,  it  seems  reasonable  to  treat  that  quantity  as  an  objective. 
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In  addition,  usually  several  competing  and  noncommensurable  objectives  ap¬ 
pear  in  a  real-life  application,  and  thus  the  designer  is  faced  with  a  decision 
making  problem  where  the  task  is  to  find  the  best  compromise  solution  between 
the  conflicting  requirements.  Many  methods  have  been  developed  to  help  the 
engineer  make  the  right  decision  in  conflicting  situations. 

In  the  majority  of  the  publications  dealing  with  structural  optimization  the 
objective,  which  should  be  minimized,  has  been  the  structural  weight.  In  many 
civil  and  mechanical  engineering  applications,  the  manufacturing  cost  may  be 
conflicting  with  the  weight  which  does  not  reflect  the  total  cost  of  the  structure. 
Accordingly,  &om  the  purely  economic  viewpoint  the  weight  represents  only  the 
material  cost  and  therefore  additional  criteria  are  needed.  When  formulating  the 
mathematical  optimization  problem,  the  designer  must  consider  which  quantities 
are  suitable  for  measuring  the  economy  and  the  performance  of  the  structure. 
Any  quantity  which  has  a  tendency  to  improve  or  deteriorate  is  actually  an 
objective  function  in  nature.  Those  quantities  which  shoiild  only  satisfy  some 
imposed  requirements  are  not  objectives  and  can  be  treated  as  constraints.  In 
a  strict  sense,  almost  all  design  quantities  have  a  criterion  nature  rather  than 
a  constraint  nature,  because  in  the  designer’s  mind  they  usually  have  better  or 
worse  values.  As  an  example  of  a  strict  constraint,  the  structural  analysis  equa¬ 
tions  (displacement,  stresses  etc.)  can  be  mentioned.  They  represent  equality 
constraints,  whereas  regulations,  space  limitations,  available  materials  and  man¬ 
ufacturing  techniques  impose  inequality  constraints.  Tl  us  the  designer  is  usually 
faced  with  several  conflicting  and  non-commcnsurable  objectives  in  formulating 
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the  optimization  problem.  In  order  to  reduce  the  computational  and  decision¬ 
making  effort  it  is  desirable  to  decrease  the  number  of  the  objectives  as  far  as 
possible  by  moving  them  into  constraints. 

A  general  multiobjective  optimization  problem  is  to  find  the  vector  of  de¬ 
sign  variables  X  =  (®i,X2,  v?rhich  minimizes  a  vector  objective  function 

over  the  feasible  design  space  X.  It  is  the 
determination  of  a  set  of  non-dominated  solutions  (or  Pareto  optimum  solutions, 
noninferior  solutions  or  efficient  solutions)  that  achieves  a  compromise  among 
several  different,  usually  conflicting,  objective  functions.  The  non-donoinated 
solutions  form  a  set  of  solutions  in  vrhich  no  decrease  can  be  obtained  in  any 
of  the  objectives  without  causing  a  simultaneous  increase  in  at  least  one  of  the 
objectives.  This  definition  can  be  understood  graphically.  An  arbitrary  collec¬ 
tion  of  feasible  alternatives  for  a  two-objective  maximization  problem  is  shown 
in  Figure  1.1.  The  area  inside  of  the  shape  and  its  boundaries  are  feasible.  The 
axes  of  this  graph  are  the  objectives  Fj  and  F2.  It  can  be  clearly  seen  from  the 
graph  that  the  noninferior  solutions  are  found  in  the  portion  of  the  boundary 
between  points  A  and  B.  Thus,  here  arises  the  decision  making  problem  &om 
which  a  partial  or  complete  ordering  of  the  set  of  non-dominated  objectives  is 
accomplished  by  considering  the  preferences  of  the  decision  maker.  Most  of  the 
multiobjective  optimization  techniques  are  based  on  how  to  elicit  the  preferences 
and  determine  the  best  compromise  solution. 
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1.1  Literature  Review 


A  variety  of  techniques  and  applications  of  multiobjective  optimization  have 
been  developed  in  the  past  few  years.  The  earliest  work  reporting  the  consid¬ 
eration  of  multiple  objectives  in  mathematical  programming  appears  to  be  that 
of  Kuhn  and  Tucker  (1951).  Progress  in  the  field  of  multicriteria  optimization 
was  suiomarized  by  Hwang  and  Masud  (1979)  and  later  by  Stadler  (1984,1987). 
Stadler  inferred  from  his  survey  that  if  one  has  decided  that  an  optimal  design  is 
to  be  based  on  the  consideration  of  several  criteria,  then  the  multicriteria  theory 
(Pareto  theory)  provides  the  necessary  framework.  In  addition,  if  the  minimiza¬ 
tion  or  maximization  is  the  objective  for  each  criterion,  then  an  optimal  solution 
should  be  a  member  of  the  corresponding  Pareto  set,  for  only  then  does  any  fur¬ 
ther  improvement  in  one  criterion  require  a  clear  trade-off  with  at  least  one  other 
criterion.  He  also  observed  that  in  mechanics  the  most  widely  used  technique 
was  goal  programming.  Bartel  and  Marks  (1974)  developed  a  method  for  ana¬ 
lyzing  optimum  design  problems  with  multiple,  competing  objective  functions. 
The  technique  provides  a  method  for  displa3ring  the  trade-off  between  two  com¬ 
peting  design  objectives  over  the  whole  range  of  the  design  space.  The  resulting 
trade-off  curve  gives  the  designer  the  best  possible  choices  for  compromise  de¬ 
signs.  The  usefulness  of  this  method  was  demonstrated  by  applying  it  to  the 
optimum  design  of  hydrodynamic  journal  bearings.  The  competing  objectives 
were  to  minimize  the  oil  film  temperature  rise  and  to  minimize  the  oil  flow  rate 
for  the  constant  loading  case.  Duckstein  (1984)  listed  about  13  possible  sample 
multiobjective  techniques  to  show  that  a  fairly  broad  choice  of  approach  is  to 
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stnictuial  design  is  possible.  Radford,  .et.al.  (1985)  in  their  study  have  explored 
the  role  of  Pareto  optimization  in  CAD.  They  used  the  weighting,  Noninferior  Set 
Estimation  Method  (NISE),  and  the  constraint  method  for  generating  the  Pareto 
optimal,  and  suggested  some  approaches  using  the  emerging  field  of  knowledge 
based  engineering  to  formalize  the  control  and  derivation  of  meaning  from  the 
Pareto  sets. 

Rao  (1984,1987,1988)  has  done  a  considerable  amount  of  study  on  multiob¬ 
jective  optimization.  In  most  of  his  research  work  goal  prograzoming  and  game 
theory  techniques  were  used.  Rao  (1984)  studied  the  application  of  multiobjec¬ 
tive  optimization  to  structural  design  problems  involving  uncertain  parameters 
and  random  processes.  The  design  of  a  cantilever  beam  with  a  tip  mass  subjected 
to  a  stochastic  base  excitation  was  considered  for  illustration.  The  objectives  in 
the  problem  were  to  minimize  structural  mass,  maximize  the  fimdamental  natu¬ 
ral  frequency  and  minimize  fatigue  damage  time.  He  obtained  numerical  results 
using  a  number  of  methods  and  concluded  that  the  game  theory  approach  was 
superior  in  finding  a  better  optimum  solution  with  proper  balance  of  the  various 
objective  functions,  even  though  the  computation  time  was  high  compared  to  the 
other  methods  discussed  in  his  work.  It  was  also  observed  that  the  optimum  so¬ 
lutions  given  by  the  various  multicriteria  optimization  methods  can  be  different. 
Rao  (1987)  proposed  a  computational  procedure  for  solving  a  general  multiobjec¬ 
tive  optimization  problem  using  cooperative  game  theory.  He  presented  graphical 
interpretations  of  the  non-cooperative  and  cooperative  game  theory  approaches 
for  a  two-criteria  problem.  Later  he  and  his  co-workers  (1988)  formulated  the 
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design  of  actively  controlled  structures  subject  to  restrictions  on  the  damping 
parameters  of  the  closed-loop  system  as  a  multiobjective  optimization  problem 
and  solved  it  using  goal  programming  techniques.  The  purpose  of  the  control 
was  to  effectively  suppress  structural  vibrations  due  to  initial  excitation.  The 
structural  weight  and  the  controlled  system  energy  were  considered  as  the  ob¬ 
jective  functions  for  minimization.  Rao,  et.al.  (1988)  formulated  a  combined 
structural/control  optimization  problem  and  solved  it  as  a  multiobjective  opti¬ 
mization  problem.  The  authors  used  the  game  theory  approach  for  generating 
the  Pareto  optimal.  The  computational  procedure  was  demonstrated  through 
the  design  of  two  actively  controlled  truss  structures.  The  two  objective  func¬ 
tions  used  for  the  optimization  were  to  minimize  the  weight  of  the  structure  as 
well  as  the  quadratic  performance  index.  Rao  (1988)  considered  the  design  of 
vibration  isolation  systems  using  multicriteria  optimization  techniques.  The  two 
performance  indices  taken  as  the  objective  functions  were  the  integrated  values 
of  the  force  transmitted  to  the  main  mass  and  the  relative  displacement  between 
the  iTutiTi  mass  and  the  base.  He  demonstrated  the  applicability  of  global  crite¬ 
rion,  bounded  objective,  utility  function,  and  game  theory  methods  in  the  design 
of  a  multicriterion  vibration  isolation  system. 

The  multicriterion  optimization  of  elastic  stress  limited  isostatic  trusses  was 
considered,  and  a  numerical  method  (weighting  method,  constraint  method,  etc.) 
for  determining  the  Pareto  optimal  set  of  the  problem  was  developed  by  Kosld 
(1982,1983).  He  discussed  the  basic  principles  of  multicriterion  optimization 
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and  presented  Pareto  optima  using  the  weighting  as  well  as  the  minimax  ap¬ 
proaches  (1981,1984).  Hyperstatic  trusses  and  plane  frames  were  designed  where 
the  weight  and  some  chosen  nodal  displacements  of  the  structure  were  taken 
as  design  criteria.  He  also  suggested  that  generating  the  whole  Pareto  set  was 
not  economical,  therefore  to  reduce  the  number  of  Pareto  optima  generated,  one 
could  develop  interactive  algorithms  where  only  a  few  Pareto  optimal  solutions 
were  generated.  Later  in  (1984)  he  developed  an  interactive  design  method  for 
truss  optimization.  Koski  (1987)  proposed  the  norm  methods  based  on  the  scalar- 
ization  of  the  original  multicriterion  problem  by  using  the  Ip-norm  method  also 
called  the  global  criterion  method.  In  addition,  an  alternative  approach  which, 
instead  of  scalarization,  reduces  the  dimension  of  the  multicriterion  problem  was 
proposed.  He  called  this  method  the  partial  weighting  method.  Throughout  the 
article  several  illustrative  truss  examples  were  presented.  In  addition  to  interac¬ 
tive  multiobjective  optimization  methods  Diaz  (1987)  proposed  sensitivity  anal¬ 
ysis  as  a  useful  tool  for  information  regarding  tradc-ofis  in  the  objectives.  Hajela 
and  Shih  (1990)  proposed  a  ininimiim  variant  of  the  global  criterion  approach  to 
obtain  solutions  to  multiobjective  optimum  design  problems  involving  a  mix  of 
continuous,  discrete,  and  integer  design  variables.  The  method  combined  a  dis¬ 
crete  variable  variant  of  the  global  criterion  approach  with  a  branch-and-bound 
strategy. 

El-Sayed,  et.al.  (1989)  demonstrated  the  formulation  of  the  nonlinear  goal 
optimization  problem  and  demonstrated  the  use  of  the  method,  as  a  design  tool. 
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Sandgren  (1989)  introduced  a  means  by  which  uncertainty  in  the  optimal  de¬ 
sign  of  structures  could  be  considered.  A  goal  programming  formulation  was 
generated  for  a  three-bar  truss  problem  wherein  uncertainty  in  both  the  load 
magnitude  and  direction  was  considered.  He  demonstrated  that  nonlinear  goal 
programming  could  be  a  valuable  tool  for  designing  under  uncertainty  conditions. 
Also  this  approach  significantly  reduces  the  complexity  of  the  solution  process. 
Dhingra,  et.al.  (1989)  presented  fuzzy  optimization  techniques  and  summarized 
their  work  as  an  useful  tool  during  the  initial  stages  of  the  conceptual  design  of 
engineering  systems  where  the  design  goals  and  the  constraints  have  not  been 
clearly  defined.  Saravanos  and  Chamis  (1990)  developed  a  methodology  for  light 
weight  low  cost  composite  structures  of  improved  dynamic  performance.  Pre¬ 
liminary  applications  on  a  cantilever  composite  beam  illustrated  that  only  the 
proposed  multiobjective  optimization  improved  all  objectives  simultaneously  as 
opposed  to  single  objective  functions.  Tseng  and  Lu  (1990)  proposed  a  minimax 
multiobjective  optimization  model  for  structural  optimization.  They  used  the 
three  typical  multiobjective  optimization  techniques,  goal  programming,  com¬ 
promise  programming  and  the  surrogate  worth  trade-off  method,  to  solve  truss 
problems.  The  main  purpose  of  their  work  was  to  apply  the  multiobjective  op¬ 
timization  techniques  to  the  selection  of  system  parameters  and  to  solve  large 
scale  structural  design  optimization  problems. 

1.2  Types  of  Multiobjective  Techniques 

Pareto  optimality  serves  as  the  basic  multicriteria  optimization  concept  in 
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virtually  all  of  the  previous  literature.  Pareto  presented  a  qualitative  definition  of 
the  optimality  concept  for  economic  problems  with  several  conflicting  objectives. 
The  Pareto  optimal  is  stated  in  simple  words  as,  “A  vector  X*  is  Pareto  optimal 
if  there  exists  no  feasible  vector  X  which  would  decrease  some  objective  function 
without  causing  a  simultaneous  increase  in  at  least  one  objective  function”.  These 
solutions  are  also  called  noninferior  solutions  or  non-dominated  solutions. 

Nearly  all  of  the  definitions  of  optimality  in  use  for  multiobjective  optimiza¬ 
tion  involve  some  sort  of  scalarization  of  the  vector  optimization  problem.  That 
is,  the  replacement  of  the  vector  problem  by  some  equivalent  scalar  minimization 
problem.  Also  the  Pareto  set  is  generally  infinite.  An  additional  use  of  scalariza¬ 
tion  is  the  selection  of  a  unique  member  of  the  Pareto  set  as  the  optimum  for  the 
vector  optimization  problem.  Usually,  a  problem  is  scalarized  either  by  defining 
an  additional  super  criterion  function  or  by  considering  the  criteria  sequentially. 

The  basic  three  techniques  used  for  generating  noninferior  solutions  are:  the 
weighting  method,  the  Noninferior  Set  Estimation  method  (NISE)  and  the  con¬ 
straint  method.  Balachandran  and  Gero  (1987)  have  discussed  the  relative  merits 
and  demerits  of  these  three  methods.  These  methods  also  come  tmder  the  cate¬ 
gory  of  nonpreference  techniques.  Although,  one  should  note  that  some  nonpref¬ 
erence  techniques  use  preference  techniques  concepts  repeated  over  a  number  of 
different  parametric  values  to  generate  the  entire  Pareto  set.  Figure  1.2.1  shows 
the  procedure  used  in  multiobjective  optimization  in  structural  optimization.  In 
the  following  sections  these  methods  are  briefly  disctrssed  with  goal  programming, 
game  theory,  global  criterion  etc. 
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Weighting  Method 


This  technique  is  based  on  the  preference  techniques  of  prior  assessment  of 
the  waghts.  It  transforms  the  multicriteria  function  to  a  single  criterion  function 
through  a  parameterization  of  the  relative  weighting  of  the  criteria.  The  problem 
can  be  stated  as: 


Minimize 


subject  to 


F(X)  =  2  WiFt(X) 


k=l 


gi{X)>0 


It  is  usually  assumed  that 


(1) 

(l.o) 


0  <  Wfc  <  1,  S  W*  =  1  (1.6) 

where  X  are  the  design  variables,  and  W  are  the  weights  for  each  of  the  q  ob¬ 
jective  functions.  With  the  variation  of  the  wdghts  the  entire  Pareto  set  can  be 
generated.  Since  the  restilts  of  solving  an  optimization  problem  can  vary  signifi¬ 
cantly  as  the  weighting  coefficients  change,  and  since  very  little  is  usually  known 
about  how  to  choose  these  coefficients,  a  necessary  approach  is  to  solve  the  same 
problem  for  many  dififerent  values  of  weighting  factors.  However,  since  the  shape 
and  distribution  characteristics  of  the  Pareto  set  are  unknowns,  it  is  impossible 
to  determine  beforehand  the  nature  of  the  variations  required  in  the  weights  so 
as  to  produce  a  new  solution  at  each  pass.  I'he  second  important  disadvantage 
of  the  method  is  that  it  will  not  identify  Pareto  solutions  in  a  nonconvex  part  of 
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the  set.  Still  conhonted  with  these  Pareto  solutions,  the  designer  must  choose 
among  them  by  some  other  means. 

Noninferior  Set  Estimation  (NISE)  Method 

The  NISE  method  (Cohon,  1978)  is  an  extension  of  the  weighting  method 
with  a  mechanism  for  quickly  converging  onto  the  Pareto  set.  In  addition,  the 
accuracy  of  the  approximation  can  be  controlled  in  the  NISE  method  through 
a  predetermined  error  criterion,  which  is  compared  to  the  maximum  possible 
error  at  every  iteration  of  the  method.  The  NISE  method  operates  by  finding  a 
number  of  noninferior  extreme  points  and  evaluating  the  properties  of  the  line 
segments  between  them.  For  example  if  two  noninferior  extreme  points  have 
been  found,  then  the  segment  between  them  is  feasible  and  it  may  or  may  not 
be  inferior.  If  the  line  segment  is  noninferior,  then  moving  in  a  direction  out 
from  the  line  segment  is  infeasible.  If  the  line  segment  is  inferior,  then  there  are 
noninferior  points  in  the  outward  direction.  Noninferior  points  are  found  in  the 
NISE  method  through  the  use  of  the  weighting  method. 

Constraint  Method 

The  constraint  method  is  discussed  in  more  detail,  since  Pareto  solutions 
are  generated  using  this  technique  for  a  few  examples  in  the  report  for  compar¬ 
isons  with  the  proposed  generalized  compound  scaling  techmque.  The  constraint 
method  is  a  technique  which  transforms  a  multicntena  objective  function  into  a 
single  criterion  by  retaining  one  selected  objective  as  the  primary  criterion  to  be 
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optimized  and  treating  the  rem£uning  criteria  as  constraints.  The  multicriteria 
problem  is  thus  stated  as: 


Minimize 


Fi{X) 

(2) 

subject  to  the  constraints 

^  *  =  2, 3,  ...,g 

(2..) 

>  0  »  =  l,2,...7n 

(2.4) 

and  the  side  bounds  on  the  design  variables  as 

«i  >  »  =  l,2,...,n 

(2.C) 

where  are  parametrically  varied  target  levels  of  the  q  —  1  objective  functions. 
Each  vector,  or  constraint  set,  B  of  the  various  fr,  will  produce  one  Pareto  solution. 
As  in  the  weighting  method,  many  different  combinations  of  values  for  each  6, 
must  be  examined  in  turn  to  generate  the  entire  Pareto  set.  The  constraint 
method  provides  direct  control  of  the  generation  of  members  of  the  Pareto  set 
and  generally  provides  an  efficient  method  for  defining  the  shape  of  the  Pareto 
set. 


Multiobjective  optimization  that  incorporates  preferences  is  described  below. 
The  commonly  used  methods  imder  this  are:  goal  programming,  game  theory, 
global  criterion,  and  the  surrogate  worth  trade-off  method. 
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Goal  Programming 


In  goal  programming  the  decision  maker  is  required  to  specify  goals  for  each 
objective  that  he  wishes  to  attain.  Goals  are  the  quantitative  values,  considered 
as  additional  constraints  in  which  new  variables  are  added  to  represent  deviations 
from  the  predetermined  targets.  The  objective  function  specifies  the  deviations 
&om  these  goals  and  priorities  for  the  achievement  of  each  goal,  in  quantitative 
terms.  Thus,  the  goal  programming  formulation  of  the  multiobjective  optimiza¬ 
tion  problem  leads  to: 


Minimize 

+  P>i 

subject  to 


(3) 


where  hj  is  the  goal  set  by  the  designer  for  the  /th  objective,  and  dj'^  and  dj  are, 
respectively,  the  under-  and  over-  achievement  of  the  yth  goal.  The  value  of  p  is 
based  upon  the  utility  function  chosen  by  the  designer.  Generally,  p  =  2  which 
is  an  Euclidean  metric  value.  The  disadvantages  of  this  method  are  the  selection 
of  “correct”  or  “valid”  levels  of  d"^  and  d~,  i.c,  which  requires  knowledge  of  the 
individual  tniTiiTnit  of  the  objective  function  which  is  not  easy  to  achieve  with 


nonconvex  problems. 
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Game  Theory 


In  game  theory,  the  multiobjective  optimization  problem  is  viewed  as  a  game 
problem  involving  several  players,  each  corresponding  to  one  of  the  objectives. 
The  system  is  considered  to  be  tmder  the  control  of  these  intelligent  adversaries, 
each  seeking  to  optimize  his  own  gain  at  the  expense  of  his  opponents,  using  all 
the  available  information.  In  this  approach,  starting  with  Xq  as  a.  starting  point, 
q  single  criterion  optimization  problems  are  solved  to  obtain  optimum  solutions 
Xi*,i  =  1,2,  ...,9.  With  the  help  of  Xi*  the  elements  of  the  matrix  [P]  are 
evaluated  as 

FiiXi*)  F2(Xi*)  ...  F,{Xi*y 

Fi{X2*)  FiiXi*)  ...  F^iXi*) 

[P]  =  ’  (4.a) 

« 

F,(V)  fjW)  n(V). 

Fkom  this  matrix  the  largest  element  of  the  t  th  column,  /'««  is  selected  from  each 
coliunn  (worst  value).  Then  a  supercriterion  or  bargaining  model  S  is  constructed 

as 

s  =  U^i^i{Fin-  Fi{Xc*)]  (4.6) 

where  Xc*  represents  the  solution  (Pareto  optimal  solution)  of  the  following  prob¬ 
lem. 


Mimmize 


F.{c,X)  =  c.fi{X)  +  (1  -  E  Fi(X)) 

i=l  i=l 
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subject  to 


gi>0,  t  =  l,2,...,m 

(4.d) 

c,-  >  0,  t  =  1,2,...,?  -  1 

(4.e) 

(4./) 

t=i 

The  supercriterion  S  defined  in  Eq  (4.b)  is  maximized, the  optimal  convex  com¬ 
bination  of  the  objective  functions  is  found,  that  is,  c*  and  the  corresponding 
optimum  solution  of  the  problem  X*  =  Xe*  is  determined. 

Surrogate  Worth  Trade-off  Method 

The  surrogate  worth  trade-off  (SWT)  method  is  an  interactive  method,  which 
can  be  applied  when  the  design  variables  are  continuous  and  the  objective  func¬ 
tions  and  constraints  are  twice  differentiable.  It  assumes  that  the  preference  of 
the  decision  maker  is  monotonic  and  can  be  exhibited  as  an  implicit  multiat¬ 
tribute  utility  function  as  the  derision  maker  systematically  compares  two  objec¬ 
tives  at  a  time.  The  SWT  method  consists  of  three  main  steps.  First,  generating 
a  representative  subset  of  non-dominated  solutions  called  promising  solutions  and 
obtaining  the  relevant  trade-off  information  for  each  generated  solution.  Second, 
interacting  with  the  derision  maker  to  elect  the  preference  structure  and  finally, 
selecting  the  best-compromise  solution  firom  the  information  obtained.  Using 
the  constraint  method  the  promising  solutions  and  also  the  associated  trade-off 
information  are  generated.  By  \ising  the  concept  duality  theory,  the  Lagrange 
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multipliers  as  well  as  the  Kuhn-Tucker  conditions,  the  trade-ofF  functions  Tii  are 
defined  as 


T. 

5/,(X) 


(5.a) 


or 


(5.6) 


Kij  is  the  generalized  Lagrange  multiplier  associated  with  the  tth  constraint. 
The  trade-off  Ay  represents  the  marginal  rate  of  substitution  between  the  tth 
and  /th  objectives.  Once  the  trade-offs  have  been  identified,  the  decision  maker 
is  supplied  with  the  trade-off  information  to  express  his  ordinal  preference. 


Global  Criterion  Method  (min-max  formulation) 


In  this  formulation  the  optimum  solution  X*  is  fotmd  by  minimizing  a  pre¬ 
selected  global  criterion  F{X)  such  as,  the  sum  of  the  squares  of  the  relative 
deviations  of  the  individual  objective  functions,  from  the  feasible  ideal  solutions. 
Thus  X*  is  found  by  minimizing 


r(x\  ^  -  Fi{X)  ^ 


(6.0) 


subject  to 

^i(A’)>0,  t=l,2,...,m  (6.6) 


where  p  is  generally  taken  as  2,  Xi*,i  =  1,2,  ...,9  is  the  feasible  ideal  solution  of 
the  i  th  objective,  and  it  is  obtained  by  minimizing  individual  criterion. 
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1.3  Scope  of  the  Work 


As  shown  from  the  description  of  the  methods,  there  are  various  drawbacks 
in  using  these  approaches  for  solving  large  scale  multidisciplinary  optimization 
problems  with  finite  element  based  analyses.  For  example,  many  single  objective 
optimization  problems  have  to  be  solved  for  generating  the  Pareto  set  in  the 
constraint  method.  Diffeient  weights  have  to  be  associated  with  the  objectives 
to  obtain  the  Pareto  set  in  the  weighting  method.  In  goal  programming,  game 
theory,  and  global  criterion  methods  individual  minima’s  of  the  single  objective 
function  are  required  which  are  difficult  to  obtain  for  a  nonconvex  problem. 
Hence,  it  can  be  seen  that  multiobjective  optimization  can  become  computer 
intensive  while  generating  the  Pareto  optimal  set. 

This  work  proposes  an  algorithm  ideally  stiited  for  multicriteria  optimiza¬ 
tion.  This  algorithm  is  based  on  the  generalized  compound  scaling  (Venkayya, 
1989-1991)  method  where  the  objectives  are  treated  as  similar  to  constraints.  A 
partial  Pareto  optimal  set  is  generated  by  the  algorithm.  The  computational  cost 
involved  in  this  method  is  much  smaller  than  the  other  methods  discussed.  This 
is  because  in  one  optimization  problem,  possible  Pareto  optimums  are  generated. 
The  knowledge  of  individual  minima’s  is  not  necessary  for  this  algorithm,  thus 
opening  a  wide  range  of  problems  that  can  be  solved  using  this  technique.  The 
algorithm  works  well,  because  the  optimal  solutions  lie  on  the  active  constraint 
boundzuy  for  structural  optimization  problems  whether  the  optimization  is  single 
objective  or  multiobjective.  In  this  work  different  objective  functions  (linear  and 
nonlinear)  were  used  to  show  the  robustness  of  the  algorithm. 
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This  report  is  divided  into  many  sections.  The  main  sections  being  the  finite 
element  formulation  for  the  plate  element,  the  sensitivity  calculations  for  the  ob¬ 
jective  functions  and  the  constraints,  the  multiobjective  optimization  algorithm 
with  mathematical  details,  and  numerical  results  with  relevant  discussions. 
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Fig.  1.2.1  Classification  of  mufticriterion  structural  design  process 
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SECTION  2 


FINITE  ELEMENT  FORMULATION 


This  section  discusses  the  finite  element  approximations  used  in  the  present  study. 
The  choice  of  elements  is  dictated  by  the  ne«d  to  perform  accurate  solutions  for  both 
thin  and  thick  plates.  The  finite  element  used  for  the  plate  bending  elements  is  a  four 
noded  ,  bilinear  displacement  element  based  on  the  Mindlin  plate  theory  (Mindlin 
1951).  A  reduced  integration  scheme  is  used  to  evaluate  the  element  matrices. 
This  results  in  a  rank-deficient  element  ,  and  must  be  stabilized  to  achieve  reliable 
behavior. 

The  theory  of  single-poim  quadrature  was  introduced  by  Hughes,  et.al.  (1978) 
and  designated  as  Ul.  This  element  was  attractive  because  of  its  simplicity,  com¬ 
putational  efficiency,  and  high  accuracy.  However  the  Ul  element  is  rank  deficient 
as  this  element  cannot  capture  the  bilinear  contributions  to  the  displacement  fields 
due  to  the  single-point  integration.  Therefore,  this  element  leads  to  some  spurious 
oscillatory  displacements  with  very  little  strain  energy  assodated  with  it.  These 
spurious  singular  modes  in  structural  mechanics  are  called  hourglass  modes  or  zero 
energy  modes.  Because  of  this  a  structure  that  appears  adequately  constrained  may 
yet  have  an  instability  that  makes  the  [K]  singular. 
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The  control  method  used  to  eliminate  instability  is  by  providing  restraint  with¬ 
out  stiffening  the  element’s  response  to  modes  of  behavior  that  are  working  well. 
The  method  adds  “hourglass  sti&ess”  to  an  element  integrated  by  one  -  point 
quadrature  .  For  a  four  noded  bilinear  rectangular  plate  element  eight  independent 
displacement  modes  can  be  identified  as  shown  in  the  Figure  2.1.  The  first  three 
are  rigid  body  modes  for  which  the  strain  energy  is  zero.  The  next  three  modes 
are  constant  strain  energy  modes.  Modes  7  and  8  are  bending  modes.  An  order 
of  one  point  quadrature  does  not  sense  these  two  modes,  as  the  Gauss  point  is  at 
the  center  of  the  element  and  hence  e*  =  e,  =  7*^  =  0  at  the  center.  Accordingly, 
the  strain  energy  for  these  two  modes  is  zero  or  almost  equal  to  zero,  hence  the 
element  displays  two  Tn«»r>iA.Tii«Trn»  (Cook,  et.o/.,  1989).  Therefore  to  provide  mode  7 
and  mode  8  with  stiffness  they  lack  under  single-point  quadrature,  a  “stabilization 
matrix  ”  is  formed  for  each  mode.  These  matrices  contain  constants  which  can  be 
chosen  such  that  a  rectangular  element  displays  the  exact  strain  energy  in  states 
of  pure  bending  (Belytschko  and  Tsay,  1983).  These  correction  stiffness  matrices 
are  added  to  the  element  stifllhess  matrix.  For  stress  computations  (especially  shear 
stresses)  the  one  point  quadrature  technique  is  the  most  accurate. 

The  bilinear  Mindlin  plate  (BMP)  elements  used  in  PROTEC  (Brockman,  ct.aZ., 
1989)  are  based  upon  uniform  reduced  integration  with  rank  corrections.  The  BMP 
bending  element  uses  single  point  quadrature.  The  rank  correction  scheme  uses  five 
generalized  hourglass  strains  to  complement  the  proper  null  space  of  the  reduced 
integrated  Mindlin  plate  gradient  operator. 
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2.1  Plate  Element 


In  the  Mindlin  plate  theory  transverse  shear  deformation  are  also  taken  into 
consideration.  Therefore  any  point  not  on  the  midsurface  is  not  governed  by  slopes 
ti;,x  zind  w,y  but  is  governed  by  rotations  and  By  of  lines  that  were  normal  to  the 
midsurface  of  the  undeformed  plate. 

u  =  -zBt  c*  =  -ztf*,*  7*v  =  +  ^v,x) 

V  —  ~~zBy  €y  —  ~^zBy^y  ^yz  “ 

7x*  =  «>,X  -  ®x 


The  state  of  deformation  is  described  by  generalized  strains, 


(7) 

(8) 


and  the  stress  state  by  generalized  forces, 

(T*  =  \NxyNy,NzyyMzjMyyMxyiQ*Z^Qy^  (®) 

In  the  stress-strain  relationships  for  the  in-surface  strains  and  curvatures,  plane 
stress  assumptions  are  used.  For  the  bilinear  plate  element  the  shape  functions  are, 

Ni  =  ;^-(s  +  + 1717  +  h»70  (10) 

Expanding  the  shape  ftmction  one  can  write 

4 
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(10.6) 


=  j(i  +  iKi  -  {) 
ff>  =  j(l  +  7)(1  +  0 
=  j(l  - 1)(1 + {) 

In  Eq  (10)  rigid  body  modes  are  denoted  by 

•‘  =  [1,1,1.11 

the  nodal  coordinates  in  the  reference  planes  are, 

{‘  =  i-i.i,i,-ij 


(10.C) 

(10.<1) 

(llo) 

(11.4) 
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Bending  strains 
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The  strain-clisplacement  matrix  evaluated  at  the  center  can  be  written  as  (Brock¬ 
man,  et.oZ.,  1989) 

/6i^  0  0 

0  0 
6a^  h\^  0 

0  0  0 

0  0  0 

0  0  0 

0  0 

^0  0  52^ 

In  Eq  (19) 

bt'*'  =  |[(ya  -  yi),  (y*  -  yi).  (y*  -  ya)>  (yi  -  y»)l  (20-«) 

ba'*'  =  |[(X4  -  Xa),  (x,  -  x,),  (x,  -  X4),  (x,  -  x^)]  (20.6) 

^  =  q  =s  0  (at  the  center  of  the  element). 


61  =  [ 


(21.a) 


II 

(21.6) 

The  element  displacement  vector 

U*  -  [u^VjV}j9xt0yy0x] 

(22) 

The  element  stiffness  matrix  is 

[K.]^llBflD][B]dv 

(23) 

K  =  B^OBA|J| 

(24) 
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D  is  the  elastic  matrix  of  the  element  which  consists  of  Dm,Df„  and  A  is  the 
area  of  the  element  and  J  is  the  Jacobian.  Theife  matrix  is  a  (  24  x  24  )  matrix. 


2.2  Stabilisation  of  the  Stiffness  Matrix 

The  plate  element  consists  of  twenty  four  degrees  of  freedom.  Ignoring  the  dz 
degree  of  freedom  the  plate  element  consists  of  six  rigid  body  modes  ,  eight  umform 
strain  modes  captured  by  the  one  point  integration  and  six  spurious  zero  energy 
modes.  The  six  zero  energy  modes  have  to  be  stabilized. 
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The  six  rigid  body  modes  are  translations  and  rotations  along  the  x,  y,and  z 
axes.  The  hour  glass  modes  can  be  represented  by  eqn(ll.d).  The  six  spurious 
modes  consist  of  five  hourglass  deformation  patterns 

u  =  h,v  =  h,w  =  hy6x  =  h,0y  =  h  (28) 

The  sixth  spurious  mode  cannot  exist  when  a  mesh  of  two  or  more  exist. 

To  control  the  kinematic  modes,  generalized  strains  are  defined  which  get  acti¬ 
vated  during  the  zero  energy  modes  and  vanish  for  rigid  body  modes.  Belytschko 
and  Tsay  (1981,1983)  defined  generalized  hourglass  strains  associated  with  each 
component  of  displacement  and  rotation.  These  are  defined  as  follows. 


91  =  93  =  7^t>,  93  =  94  =  7^^,,  95  = 


where 

7  =  h  -  -f  k^ybj]  (29) 

Equation  (29)  is  valid  because  of  the  orthogonality  conditions.  The  generalized 
hourglass  stiffiiess  associated  with  the  modes  described  in  Eq  (28)  are 


/  1  \ 


/Ei\ 

1 

E2 

O.lEtA 

Ei 

—  -t  .  1 

u 

Ea 

\Eij 

t» 

13 

\TiJ 

(30) 


These  stifihesses  are  added  to  the  original  stififhess  matrix. 
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2.3  Mass  Matrix  Formulation 


In  static  problems  the  stabilization  scheme  mentioned  (Brockman,  et.al.,  1987) 
works  pi^rfectly,  giving  good  element  behaviour.  In  dynamics,  however,  low-energy 
deformation  patterns  consisting  mainly  of  hourglassing  motions  represent  likely 
modes  of  low-frequency  oscillations.  These  low  frequency  oscillations  contaminate 
the  vibration  spectrum  which  is  of  the  most  interest. 

This  can  be  shown  by  the  following  method.  Given  a  vibration  mode  shape  d  = 
$,  the  generalized  stiffness  and  mass  associated  with  the  mode  are  the  projections: 

k  =  m  =  (31) 

and  the  corresponding  frequency  of  vibration  is 

(32) 

Since  the  hourglass  modes  are  orthogonal  to  the  constant  strain  modes  and  rigid 
body  modes,  the  stiffness  associated  with  these  modes  depends  solely  on  the  hotur- 
glass  stiffness  which  is  deliberately  made  small  to  avoid  locking  problems.  But  at 
the  same  time  the  kinetic  energy  associated  with  these  unstable  modes  is  simjiar  in 
magnitude  to  that  a.ssociated  with  the  rigid  body  and  uniform  state  motions.  This 
is  because  the  orthogonality  conditions  are  no  longer  existing  while  calculating  the 
kinetic  energy  of  the  system.  Consequently  artificial  vibrations  modes  appear  in  the 
element  spectrum  ,  intermixed  with  the  lower  fundamental  frequency  modes  which 
are  commonly  of  interest. 
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Hence  during  the  mass  matrix  formulation  he  mass  matrix  has  to  be  corrected 
by  reducing  or  eliminating  the  kinetic  energy  associated  with  unstable  modes  of  the 
element,  while  leaving  the  energy  associated  with  rigid  body  motions  and  uniform 
strain  states  imchanged  (Brockman,  et.al.,  (1987)). 

For  the  mass  properties  of  the  bilinear  element 

/( 

*  p(l,z,z*)dz 
r 

and 

H  =  NNUA 

The  kinetic  energy 

r  =  ^  /  +  v’  + 1"’)  +  2Ri{uey  -  vK)  +  ^3(^1  f  (35) 

2<fA 


(33) 

(34) 


The  consUtent  mass  matzix  formulatioa  is  as  follows  from  Eqs  (33,34): 


/Riff 

0 

0 

0 

R2ff\ 

0 

RiH 

0 

-R^ff 

0 

0 

0 

Riff 

0 

0 

0 

-R^H 

0 

RiH 

0 

V  /ijff 

0 

0 

0 

Riff) 

Assuming  a  constant  Jacobian  determinant  the  above  equation  reduces  to 


H  = 


(36) 


(37) 
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Consistent  Mass  via  the  Projection  Method 


The  kinetic  energy  of  the  hourglass  modes  is  relatively  large  compared  with  the 
hourglass  stiffness  used  for  stabilization.  Hence  the  mass  matrix  is  formulated  as 
follows,  considering  the  kinetic  energy  associated  with  the  in  displacement  u  : 

Tu=  I  l-Riu^dA  (38) 

Ja  2 

We  have  to  eliminate  the  kinetic  energy  associated  with  that  part  of  the  velocity 
field  which  is  ignored  by  the  single  point  integration  of  the  element  stifihess,  and 
then,  expand  the  velocity  field  u{x,y)  about  x  =  y  =  0.  Hence 

^*1  y)  =  0)  +  xu,x(x,  y)  +  yu,y(x,  y)  +  ...  (39.o) 

or 

+  x(5i*u)  +  yih^ti)  +  rj({'r*u)  (39.6) 

Modifying  the  kinetic  energy  based  on  the  purely  linear  part  of  the  velocity  field, 
that  is,  keeping  only  the  linear  terms  &om  Eq  (39.b) 

U{(x,y)  =  7S*u  +  z{bi*u)  +  y(62*u)  (40) 

4 

gives 

Ti  =  -tt^f  /  i2i(-s  +  x6i  +  y63)(7s  +  ®6i  +  y62)  dA]u  (41) 

2  Ja  A  4 

As  with  the  stifihess  computation  ,  the  Jacobian  determinant  is  a  constant 

=  J  ^  ^(s  +  U  +  VV  +  ^iv)d‘ndi  =  I'f I*  (42) 
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K  the  center  of  the  element  is  x  =  y  =  0  then 

jJ^xdA  =  jT  N*xdA  =  |J|(s*x)  =  0  (43) 

and  hence  terms  linear  in  x  or  y  do  not  survive  the  integration. 

J^x*dA  =  x*Hx  =  Cx»  (44.o) 

y  VA  =  y*Hy  =  Cyy  (44.6) 

^xy*dA  =  x‘Hy  =  c,y  (44.c) 

Therefore  the  kinetic  energy 

Ti  =  ^R^iL^Hu  (45) 

in  which 

A 

B  =  — sa*  +  c„b,b,*  +  Cy^bab,*  +  Cxy(biba*  +  bjb/)  (46) 

by  performing  a  similar  linearization  of  all  displacements  and  rotation  components, 
H  is  replaced  by  It  is  evident  that  the  kinetic  energy  associated  with  all  the 
five  hourglass  modes  is  identically  zero.  And  hence  modes  consisting  primarily  of 
hourglass  motions  should  not  appear  as  spurious  low-energy  vibration  modes.  The 
mass  matrix  obtained  by  projection  onto  a  linewzed  velocity  field  is  positive  semi- 
defimte,  since  zero  kinetic  energy  is  associated  with  all  pure  hourglass  patterns. 
This  mass  formulation  is  valid  for  plate  bending  alone,  since  the  singular  modes 
of  the  mass  matrix  coincide  precisely  with  those  of  the  stiShess.  The  projected 
mass  matrix  3rields  the  proper  energy  for  all  eight  elementary  states  and  zero  for 
the  hourglass  modes.  The  mass  obtained  by  single  point  quadrature  leads  to  a 
proper  energy  only  for  the  translational  rigid  body  modes.  Based  on  this,  a  single 
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point  quadrature  is  used  for  inplane  motions,  and  the  projection  method  for  the 
transverse  displacements  and  rotations. 


Mppt  — 


( BiH(i) 
0 
0 
0 

Rik 


0 

R^Hil) 

0 

-R2H 

0 


0 

0 

RxH 

0 

0 


0  Rik^ 
-Rik  0 
0  0 
Rik  0 
0  Rik  J 


(47) 


where  the  element  is  sampled  at  the  element  centroid,  N  =  ^s.  Therefore, 


16 


1  1 
1  1 
1  1 
1  1 


(48.a) 

(48.6) 


The  semi  property  of  the  single-point-integrated  mass  and  the  projection  mass  ma¬ 
trix  appears  to  present  a  potential  source  of  difficulty  in  some  methods  of  eigenvalue 
extraction,  such  as  subspace  iteration.  But  the  subspace  projection  of  the  mass  will 
remain  positive  definite  unless  one  or  more  trial  vectors  correspond  precisely  to  a 
global  deformation  mode  which  is  free  of  kinetic  energy.  This  situation  is  highly 
unlikely. 
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SECTION  3 


SENSITIVITY  ANALYSIS 

This  section  describes  the  calculations  of  sensitivity  information  for  the  system  by 
the  direct  method.  The  sensitivity  parameter  of  interest  is  the  thickness  of  the  plate. 
The  optimization  of  structures  under  static  loading  typically  requires  derivatives  of 
displacement  or  stress  constraints.  Because  element  stresses  may  be  easily  written 
aa  explicit  functions  of  the  nodal  displacements,  a  general  form  of  the  constraint  is 

^(u,t)>0  (49) 

In  Eq  (49)  u  is  the  displacement  vector  and  t  is  the  design  variable.  The  derivative 
of  g  with  respect  to  the  design  variable  t  may  be  written  as 

dt  at  ^  dt 

where  Z  is  a  vector  such  that 

The  first  term  in  Eq  (50)  is  usually  zero  or  very  easy  to  obtain,  so  the  compu¬ 
tation  of  the  second  term  is  discussed  below. 
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3.1  Displacement  Constraints 


The  static  equilibrium  for  the  structure  is 


[K]U  =  P 


(51) 


where  K  is  the  global  stifihess  matrix  and  P  is  the  load  vector.  Differentiating  Eq 
(51)  with  respect  to  the  design  variable  t  we  get 

For  load  P  not  varying  with  thickness  the  above  equation  becomes 

Premultiplying  Eq  (52.b)  by  we  obtain 

zt^^zTk-'(-i^v)  (M-c) 


This  above  equation  is  solved  once  for  each  design  variable.  Let 


k'u 

dt 

(53) 

we  know 

K  =  B'^DB 

(54) 

therefore 

K*  =  B’^d'B 

(55) 

as  only  the  D  matrix  contains  t.  Therefore 

KU  =  B'^d'BU 

(56) 
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or 


where 


K'u  = 


a  =  A 


/ 


V 


Ua. 

t 

t 

< 

3^ 

3^ 

3^ 

Qm. 

t 

9i. 


\ 


/ 


(57) 


(58) 


3.2  Stress  Constraints 

Once  the  solution  for  U'  is  complete  from  £q  (52.c)  the  stress  sensitivities  can 
be  obtained  from  the  following  sequence.  The  stress  derivative  is  derived  element 
by  element. 


ff  =  DBU 


(59.a) 


^  =  d'bU  +  DBU' 


(59.6) 
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3.3  Frequency  Constraints 
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(64.5) 


duj  _  {Uj'^K'Uj)  -  uj^Uj'^M'Uj) 
dt  ~  2u)jUj'^MUj 

_  Uj'^[K>  -  utj^M']Uj 
2u;j{Uj'^MUj) 

for  the  j**  mode  of  vibration. 

The  denominator  is  a  scalar  mdtiple  of  the  generalized  mass  for  the  mode  j, 
which  can  be  calculated  at  the  system  level.  The  product  Uj^K'Uj  is  calculated 
element  by  element  and  then  assembled  globally.  The  same  is  the  case  for  the  mass 
matrix  sensitivity. 

Concentrated  mass  added  at  each  node  will  not  affect  the  sensitivity  calculation 
of  the  mass  matrix,  as  the  concentrated  mass  value  is  a  constant.  The  total  mass 
M  =  M,  +  Me  where  Af«  is  the  mass  due  to  the  structure  and  Me  is  the  mass  due 
to  the  concentrated  mass. 


Differentiating 


,  j  dM  d[M,  +  Me] 

at 

(65.a) 

a(jif.j  a[M,\ 
at  at 

(65.5) 

where 

me]  0 
at 

(65.c) 

Hence  the  concentrated 

mass  does  not  have  any  affect  on 

the  mass  sensitivity.  Mass 

M  is  denoted  as 

M  =  J^RNN'^\J\dV 

(«6) 

(64.c) 
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were  R  is  a  function  of  the  element  density  and  thickness  .  As  R  is  a  constant  for 
an  element,  it  can  be  taken  out  of  the  Integral.  The  element  mass  matrix  can  be 
differentiated  with  respect  to  the  thickness  easily. 

Typical  procedures  for  system  level  solution  are  outlined,  sissuming  that  element-| 
level  routines  are  available  for  evaluating  the  vectors  K  U  and  M  U ,  and  the  scalar 
products  U'^K  U  and  U^M  U  as  required.  For  the  eigenvalue  problem,  for  each 
design  variable  and  mode,  the  frequency  sensitivity  can  be  summed  element  by 
element. 

(67..) 

^3^3  e=l 

mj  =  Uj'^MUj  (67.6) 


(j  is  not  summed) 


SECTION  4 


OPTIMIZATION  ALGORITHM 


This  section  discusses  the  development  of  the  scaling  algorithm  and  presents  the 
algorithm  for  multiobjective  optimization  with  a  discussion  on  reliability  calcula¬ 
tions  which  are  used  for  selecting  the  best  compromise  design. 

The  algorithm  proposed  is  based  on  optimality  criteria  methods  (Venkajrya, 
1986).  The  optimality  criteria  method  has  the  potential  to  extend  optinuzation 
to  problems  with  thousands  of  variables  and  constraints.  One  of  the  significant 
elements  of  the  optimality  criteria  method  is  the  concept  of  scaling.  Scaling  im¬ 
plies  ffhanging  the  Variables  with  the  objective  of  reaching  the  constraint  boundary. 
The  design  variables  can  be  brought  to  the  constraint  boundaries  either  by  adding 
differential  quantities  or  by  multiplying  the  variables  by  scale  factors.  The  lat¬ 
ter  is  called  scaling.  This  procedure  was  originally  proposed  for  constraints  which 
could  be  scaled  in  a  single  step.  Later  it  was  generalized  for  nonlinear  functions. 
In  both  cases,  however,  a  simple  scaling  algorithm  was  used,  and  it  was  adequate 
for  structural  optimization  problems.  In  simple  scaling  a  dominant  constraint  can 
be  identified,  and  scaling  with  this  constraint  generally  brings  the  design  variables 
into  the  feasible  region.  This  is  not  necessarily  the  case  in  general  mathematical 
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optimization  problems,  and  a  compound  scaling  is  necessary  in  order  to  extend 
the  scope  ot  the  optimality  criteria  method.  Another  significant  element  of  the 
optimality  criteria  method  is  modification  of  the  variable  vector  at  the  constraint 
boundary  by  directly  invoking  the  optimality  as  defined  by  the  Kuhn- Tucker  con¬ 
ditions.  The  new  techniques  which  were  developed  for  multiobjective  optimization 
based  on  scaling  are  discussed  in  Section  4.2.  The  following  sections  explicitly 
discuss  the  development  of  the  scaling  algorithm. 

4.1  Development  of  Scaling  Algorithm 

The  first  objective  of  scaling  is  to  bring  all  the  violated  constraints  to  the  feasible 
region  or  to  the  constraint  boundary.  If  all  the  constraints  are  in  the  feasible 
region,  then  scaling  can  be  used  to  bring  some  of  the  constraints  to  the  constraint 
boundary.  The  mathematical  basis  for  scaling  can  be  derived  from  the  first  order 
approximation  of  a  Taylor’s  series  (Venkayya,  1971,  1986). 

4.1.1  Scale  Factor  Derivation  for  Simple  Scaling 

Simple  scaling  is  done  with  respect  to  the  most  prominent  constraint  (Figure 
4. 1.1.1)  The  scaling  procedure  can  be  explained  with  the  help  of  two  designs  rep¬ 
resented  by  the  two  variable  vectors  x  and  x.  The  relationship  between  the  two 
variable  vectors  is  given  by 

X  =  Ax  (68) 
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where  A  is  a  single  scalar  parameter  which  will  be  referred  to  as  a  scale  factor. 
(A  >  0).  If  dx  is  the  difference  vector  between  the  two  designs,  then  one  can 
write 


dx  =  X  —  X  =  (A  —  l)x  (69) 

Also  if  g  and  g  are  the  response  quantities  respectively  in  the  two  designs,  then 
a  change  in  response  can  be  represented  by 

dg  =  g-g  (70) 

Now  from  the  definition  of  the  total  differential  (  first  order  approximations  of 
the  Taylor’s  Series)  the  following  relationship  can  be  written 

,  dg  dg  j  dg 

dg  =  — dxi  +  — dx2  +  ...  -r  — dxn 
dXi  0X2  dXn 

Then  dg  can  also  be  written  as  (from  Eqs  69  and  71) 

^^  =  (A-l)EiVii®i  (72) 

1=1 

where  Nn  =  Dividing  Eq  72  by  g  we  have, 

-  =  (A  -  l)i;  —  (73) 

9  i=ii  9 

An  examination  of  Eq  73  presents  two  interesting  cases. 


Case  1: 


<0 

9 


In  this  case  a  new  paremieter  fi  is  defined  as 


(74) 


Then  Eq  73  can  be  written  as 


Then 


Solving  for  the  scale  factor  A 


A  = 


fijN 

+  y-jN  -  1 


(84) 


This  scaling  technique  was  originally  developed  by  Venkayya  (1971),  without  la¬ 
belling  various  terms  as  ^  and  fi  parameters.  Nevertheless,  the  concept  has  been  in 
use  for  the  last  two  decades.  Substituting  A,  and  /?  expressions,  we  get 


A- 

s-r.U;^i-9 

(85) 

A 

(86) 

(87) 

(88) 

Equation  88  is  the  reciprocal  approximation  used  in  mathematical  programming 


techniques,  which  is 


J(»)  =  S(*)  +  E  i— 


(89) 


This  is  the  form  Schmit  and  Farshi  presented  in  1974  as  the  reciprocal  constraint 
approximation,  which  is  a  very  popular  technique. 


Case  2: 


>  0 


Now  the  parameter  /ijp  is  defined  as 


^  Niixi 
f=l  ^ 


(90) 


(91) 
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Then  the  scale  factor  A  can  be  written  as 


^  ^  P  +  II jp  -  1 


HP 


Decoding  the  above  equation  (92)  we  get 


i  .  Efe  _  ^ 

A  =  J - f - 


(92) 


(93) 


Simplifying 


g  +  E  fexi  -  g 


g{x)  =  ^(a!)  +  -  1) 

=  i/(®)  +  S  “  *•) 


(94) 


(95.a) 


(95.6) 


the  above  equation  is  the  linear  approximation  used  in  mathematical  programming 
techniques  (Schmit  and  Farshi,  1974). 


An  examination  of  Eqs  84  and  92  reveals  some  interesting  facts: 


1.  In  Case  1  the  scale  factor  is  inversely  proportional  to  the  target  response  ratio, 
and  in  Case  2  it  is  directly  proportional  to  /?. 

2.  The  response  of  the  system,  </,  and  the  response,  Nn  can  be  determined  from 
an  analysis  of  the  system  for  a  given  variable  vector  z.  The  target  response  (or 
desired  response)  can  be  determined  from  the  constraint  definition.  Then  the 
target  response  ratio  P,  and  the  parameter  fi  are  known.  Then  the  scale  factor 
A  can  be  determined  explicitly  for  any  type  of  structure  and  constraints. 
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4.1.2  Compound  Scaling 


In  compound  scaling,  scaling  is  done  with  respect  to  multiple  constraints  Figure 
4.1. 2.1.  The  scale  factors  derived  from  the  first  order  approximations  of  a  Taylor’s 
series  Eqs  84  and  92  represent  an  interesting  generalization  for  nonlinear  functions. 
A  formula  involving  linear  and  nonlinear  interaction  was  written  by  Venkayya  (1989) 
in  which  it  was  stated  that,  when  the  constraint  is  nonlinear  the  scaling  formula 
can  be  represented  by 

A  a  (i)s  (96) 

where  the  parameter  m  is  a  measure  of  nonlinearity. 

The  basis  for  generalization  of  the  scaling  algorithm  can  be  established  by  ex* 
amining  the  and  parameters  of  two  simple  constraint  functions. 

=  +  ^  +  (97) 

*1  ®2  ®n 

g2  =  CiX\  +  C2Z3  +  —  +  CnZ„  <  ^  (98) 

where  *i,X2,  are  the  variables  and  ci,C29...,Cn  are  a  set  of  positive  constants. 
The  n  parameters  in  these  two  cases  are  defined  as 

=  =1  (99) 

»=i 

=  i  (i»») 

i=i  92 

The  measure  of  nonlinearity  in  both  cases  is  1,  and  it  represents  the  linear  case. 
The  parameters  0  are  defined  as 

A  =  S  (101) 

91 
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92 

It  is  easy  to  show  that  the  exact  scale  factors  in  these  two  cases  are  simply 


(102) 


(103) 


A2P  =  P2 


(104) 


The  first  condition  represents  Case  1,  and  the  second  represents  Case  2  in  the 
simple  scaling  technique  discussed  earlier.  A  similar  examination  was  made  of  two 
nonlinear  functions 

(105) 


9i  -  +  T";;:  +  —  +  — ^ 


*2"*  ®n”* 


g2=C\Xi  +  caajj"*  +  ...  +  CnXn  <  ^2 


(106) 


The  n  parameters  for  these  two  cases  are 


MiJV  =  m 

(107) 

H2P  =  m 

(108) 

The  meastire  of  nonlinearity  in  these  two  cases  in  m.  It  is  easy  to  again  show  that 
the  exact  scale  factors  for  these  two  cases  are 

Aijv  =  (109) 

Aap  =  (110) 

For  constraints  with  both  linear  and  nonlinear  terms  it  can  be  shown  that 


1  < 


(111) 
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The  scale  factors  given  in  equation  109  and  110  are  far  superior  approximations  than 
those  given  in  Eqs  84  and  92.  However,  they  are  not  exact  when  the  summation  in 
the  fi  parameter  definition  contains  both  positive  and  negative  terms.  In  such  cases 
two  scale  factors  are  defined  for  each  constraint.  These  correspond  to  Case  1  and 
Case  2  discussed  earlier.  For  example,  the  two  scale  factors  for  the  constraint 
can  be  defined  as 


■V  = 

A,7>  = 

The  definition  of  Hjff  and  fijp  in  Eq  112  and  113  is  as  follows; 

A  NijXi 
«=1  Si 


(112) 

(113) 


and 


^  NijXi 


(114) 


(115) 


where  Eq  114  contains  only  negative  terms  and  Eq  115  contains  only  positive  terms. 
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4.1.3  Constraint  Approximations 


The  conservatiye  approximation  used  the  above  concept  in  building  the  approx¬ 
imate  ftmction.  The  constraint  approximation  is  defined  as  (Starnes  and  Haftka, 
1979) 


n 


9j  —  9i  +  (Aj  ~ 

isz\ 

(116) 

where 

Gj  —  NijXi  <  0 

^jN 

(117.a) 

Gj  =  1  NijXi  >  0 

(117.6) 

Dividing  Eq  116  by  gj 

9j  i=i  9j 

(llS.tt) 

^  =  1  -h  (Ay  —  l)MjGj 

(118.5) 

If  fij  is  positive,  then  it  becomes  ftjp  and  Gj  becomes  1.  Therefore 

0=1  +  (Ayp  -  l)pjp 

(119.o) 

and 

.  0  +  HP  -  1 

\jP  — 

HP 

(119.5) 

which  is  the  same  as 

the  general  equation  derived  for  nonlinear  constraints  given 

in  Eq  113  (when  fijp 

is  equal  to  1).  If  ^y  is  negative,  then  it  becomes  fijn 

and  Gj 

becomes  -r^.  Therefore 

(120.a) 

and 

A.  - 

/?  +  HN  -  1 

which  is  the  same  as  £q  112  {fijif  is  equal  to  l). 


(120.fc) 
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4.2  Mutiobjective  Compound  Scaling  (MCS)  Algorithm 

Multiobjective  optimization  was  solved  by  modifying  the  generalized  compound 
scaling  algorithm  (Venkayya,  1991)  for  multiobjective  optimization  problems.  This 
technique  generates  a  partial  Pareto  optimal  set.  For  a  multiobjective  problem 
there  exists  several  Pareto  optima.  Since  the  determination  of  all  the  Pareto  optima 
of  a  mtiltiobjective  problem  is  computationally  expensive,  a  subset  of  the  Pareto 
optimal  points  is  sufficient  for  finding  the  solution.  As  stated  by  Rao  (1988)  and 
observed  by  many  other  researchers  the  optimum  solution  generated  by  the  various 
multicriteria  optimization  methods  can  be  different  from  each  other,  hence  a  solution 
concept  or  procedure  should  be  defined  on  the  basis  of  other  attributes  such  as  the 
mathematical  basis  of  the  method,  its  generality  and  the  quality  of  the  final  solution. 
Keeping  this  in  mind  the  generalized  compound  scaling  algorithm  for  multiobjective 
or  vector  optimization  is  built. 

The  MCS  algorithm  is  based  on  scaling  techniques.  The  basic  idea  of  any  scaling 
algorithm  being  to  derive  a  set  of  scale  factors  for  the  variables  such  that  the  design 
can  be  brought  to  the  constraint  surface  in  one  or  more  steps  from  anywhere  in  the 
n-dimensional  design  space  by  multiplying  the  current  values  of  the  variables  and 
the  scale  factors.  The  implication  being  that  the  optimum  lies  on  the  constraints 
surface  or  intersection  of  the  constraints.  To  generalize  the  scaling  techmque  the 
objective  functions  are  also  treated  as  additional  constraints,  as  both  the  constraints 
and  the  objective  function  are  functions  of  the  same  variables.  In  multiobjective 
optimization  it  is  well  known  that  the  Pareto  optimal  set  lies  on  the  intersection  of 
the  objective  functions  contours,  hence  treating  the  objective  functions  in  the  same 
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way  as  the  constraints  seems  a  logical  approach.  So  once  an  optimum  constraint 
surface  or  intersection  of  constraint  surfaces  is  reached,  the  MCS  technique  tries  to 
reduce  one  of  the  objective  functions  depending  on  the  objective  function  gradient 
information.  The  difficulty  in  treating  the  objective  functions  as  constraints  is 
that  the  constraints  have  a  target  to  satisfy  and  the  objective  functions  do  not. 
Hence  to  make  the  objective  fimctions  look  like  constraints,  pseudo  target  values 
for  the  objective  functions  are  generated  so  that  they  can  be  treated  as  additional 
constraints  which  intersect  other  constraints.  Linear  as  well  as  nonlinear  objective 
functions  were  considered  to  demonstrate  the  robustness  of  the  algorithm. 

Multicriteria  Problem 

The  problem  of  multiobjective  optimization  subject  to  multiple  constraints  can 
be  posed  as: 

Minimize 

f  =|f,(x).F,{x) . fi(x)l  (121) 

subject  to  the  constraints 

•••>*»»)  S  9i  j  ~  l,2,...,o. 

~  j  =  o  +  l,...,7n. 

and  the  side  bounds  on  the  design  variables  as 

®*  >  *  =  1»  2, ...,  n. 
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(122.a) 

(122.6) 

(123) 


where  F  are  the  conflicting  objectives  to  be  minimized,  q  is  the  total  number  of 
objectives,  n  is  the  number  of  design  variables  and  m  is  the  total  nximber  of  con¬ 
straints. 

MCS  Algorithm 

The  technique  makes  use  of  the  active  set  constraint  strategy  hence  reducing 
the  computational  efibrt  considerably  (Venkayya  and  Tischler,  1989).  The  two  most 
important  parameters  are  the  target  response  ratio  and  the  characteristic  function. 
They  are  denoted  by  /3  and  /i,  respectively.  In  general,  to  put  it  in  a  broad  sense, 
the  ^'a  are  the  scale  factors  magnitudes  for  the  new  design  vector,  and  the  n'a 
are  the  direction  vectors  for  the  magnitudes.  The  Figure  showing  the  sequence  of 
operations  is  given  in  Figure  4.2.1. 

First,  0  (target  response  ratio)  parameters  are  computed  for  all  the  constraints 
using  the  following  expression  (same  as  Eq  82). 

/?,  =  ^  »  =  l,2,...,m  (124) 

9i 

If  gi  and  arc  positive,  then  the  following  statements  are  valid.  0i  =  1  represents 
the  constraint  surface,  /?,•  >  1  represents  the  feasible  region,  and  0i  <  1  represents 
the  infeasible  region. 

After  arranging  the  constraints  with  the  0^s  in  ascending  order,  the  constraints 
with  the  lowest  0  values  being  the  most  critical  are  selected  for  the  ictive  constraint 
set  (about  10  to  15)  and  denoted  by  a.  Next,  the  characteristic  function  parameter 
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H  is  computed  for  each  of  the  selected  a  constraints  and  the  objective  functions. 
The  /i  parameter  consists  of  two  parts  based  on  the  sign  of  the  individual  terms  in 
the  summation.  This  basically  stores  the  gradient  information  of  the  constraints 
and  the  objective  function. 


For  all  negative  mj  terms 


HN  =  -  E  i  =  h  2,  a  +  g 


(125.a) 


For  all  positive  /iij  terms 


HjP  =  5^  liij  j  =  1, 2, a  +  q 
1=1 


where 


NijXi 


y  ^  Ij  2j  ^  q 


(125.6) 


(125.C) 


where  zj  is  a  vector  containing  the  values  of  a  selected  constraints  and  the  q  objective 
functions.  Nij  is  the  ordered  gradient  matrix  corresponding  to  the  selected  active 
constraint  set  a  in  the  order  of  their  0  values  which  are  sorted  in  ascending  order 
and  derivatives  of  the  objective  functions  with  respect  to  the  design  variables. 


To  make  the  objective  fimctions  look  like  constraints,  one  may  write 

Fi{x)<Ti  »  =  l,2,...,g 

where  7%  are  pseudo  targets  for  the  objective  functions.  Two  distinct  cases  are 
identified  for  assigning  the  P  values  for  the  objective  functions. 
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Case  1: 


When  the  design  is  away  from  the  constraint  surface, 


then 


1  +  a 

(126.a) 

)9jb  =  '!'  >  HP 

(126.6) 

^  1 

^  HN  <  HP 

(126.c) 

where  A  =  s  +  l,s  -1-  2,  ...q,  and  /3i  is  the  target  response  ratio  for  the  most  violated 

or  most  active  constraint  from  the  active  constraint  set  a.  The  a  is  calculated  using 

the  information  provided  from  the  corresponding  characteristic  function  ft  for  the 

same  selected  constraint. 

1 

7  =  — 

H 

(127.a) 

a  =  7  HN  >  HP 

(127.6) 

1 

a  =  -  HN  <  HP 

(127.c) 

where  m  is  the  larger  value  of  the  positive  (/*ip)  or  negative  {h\a)  for  the  selected 
constraint. 

Case  2: 

In  the  second  case  when  the  design  is  on  the  constraint  boundary,  the  ^  pa¬ 
rameter  of  the  objective  functions  are  weighed,  and  the  objective  function  which 
has  the  smallest  value  is  selected  for  scaling  down.  Then  the  objective  function 
selected  is  scaled  down  by  multiplying  or  dividing  the  current  values  of  the  variables 
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by  a  factor  depending  on  whether  the  gradient  of  the  objective  function  is  positive 
or  negative  with  respect  to  that  variable.  The  value  of  the  factor  is  generally  taken 
as  0.8. 


The  next  step  is  to  calculate  the  scale  factors.  Scale  factors  A  depend  on  the 
values  of  the  function  characteristic  parameters  and  target  response  ratios.  The 
procedure  is  as  follows: 


AyAT  = 


fijN  >  HP 


(128.a) 


Ayp  = 


HN  <  HP 


(128.6) 


The  following  conditions  are  used  in  constructing  the  scale  factors  matrix. 

(129  “) 


< 

II 

< 

IHi<0 

Aij  =  1.0 

Hi  ^  0 

Aiy  =  Ayp 

Hi  >  0 

Ay  =  1.0 

Mij  <  0 

HN  <  HP 


(129.6) 


This  results  in  s  +  9  possible  scale  factors  for  each  variable.  The  strategy  in  com¬ 
pound  scaling  is  to  select  a  mix  of  the  scale  factors  from  various  columns  of  the 
scale  factor  table  with  the  object  of  approaching  the  intersecting  points  of  the  con¬ 
straints.  The  relevant  scale  factors  are  determined  with  the  help  of  a  scale  factor 
assignment  procedure.  The  underlying  concept  is  similar  to  pivoting  for  identifying 
prominent  variables.  The  scale  factors  assignment  table  is  computed  as  follows 


Ui  =  \^^\^HN  Hi 

tij  =  0  /I 


yy<0  I 

Hi  >  0  J 


HN  ^  HP 


(130.0) 
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or 


Uj  =  0  /iti  <  0  J 


<  HP  (130.6) 


The  values  of  the  entries  in  the  scale  factor  assignment  table  vary  from  zero  to  one. 
Now  all  the  information  is  available  for  the  selection  of  an  appropriate  scale  factor 
for  each  variable  from  the  list  of  s  +  9  scale  factors.  The  appropriate  scale  factor 
for  the  variable  is  selected  based  on  the  following  principles: 


(i)  From  the  scale  factor  assignment  matrix  select  the  largest  entry  in  the  t**  row. 


For  Example:  If  Uj  is  the  largest  entry  in  the  scale  factor  assignment  matrix, 
then  the  appropriate  scale  factor  for  the  variable  is  the  A.ij  from  the  scale  factor 
matrix. 


(ii)  If  there  is  more  than  one  entry  equal  to  the  largest  entry  in  that  row,  then 
the  appropriate  scale  factor  corresponds  to  the  smallest  Pj  of  the  equal  entries. 
Now  this  process  of  selecting  the  appropriate  scale  factor  is  applied  to  all  the 
variables.  The  selected  scale  factor  for  the  variable  is  denoted  by  r<.  The 
scaled  variables  for  xj  are  computed  by 

(131) 


—  _  -.—  .old 

Zi  =  TiXi 


t=  1,2. 


The  compotmd  scaling  algorithm  simply  drives  the  design  to  the  intersection  of 
the  constraints,  if  such  an  intersection  exists.  This  can  be  illustrated  well  with  an 
example  which  consists  of  two  design  variables  and  two  constraints  for  a  bicriteria 
optimization  problem.  The  problem  statement  is  as  follows: 
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Demonstration  Problem 


Minimize 


Fi  =  4x1  + 


F2  = 


0.004 

X2  +  0.25xi 


Subject  to 


91  = 


1.6 


§2  =  0.2[^^  + 


xj  -|~  0.25x1 
2 


3xi  xj  +  0.25xi 


<  1.0 
1  <  1.0 


(132,a) 

(132.6) 


(133.a) 

(133.6) 


Start  with  an  initial  design  xo  =  (0.1, 1.0)  which  are  the  lower  boimds  of  the  design 
▼ariables.  It  can  be  seen  from  Figure  4.2.2  that  in  one  scaling  the  design  reaches 
the  intersection  of  the  two  constraints  gi  and  g2.  This  comes  about  because  of 
the  resultant  vector  direction  of  the  gradients  of  the  two  constraints,  and  the  scale 
factor  is  calculated  as  discussed  above.  At  point  A  in  Figure  4.2.2  it  can  be  seen 
that  moving  in  the  direction  of  objective  fimction  Fi  for  miniimzation  will  violate 
the  constraints.  Hence  at  point  A  the  designs  are  scaled  by  a  factor  0.8  with  respect 
to  F2  depending  on  the  gradients  of  the  second  objective  function.  This  is  point  B 
in  Figure  4.2.3.  From  point  B  the  design  moves  to  the  intersection  of  gi^F\yF2  to 
point  C  shown  in  Figure  4.2.3.  The  dark  highlighted  curve  shows  the  minimal  curve 
(  Pareto  optimum  ).  In  this  procedure  the  design  moves  into  the  feasible  region  and 
then  back  to  the  constraint  boundary. 
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4.3  Reliability  Criteria 


In  multiobjective  optimization  there  are  a  number  of  Pareto  optimal  design, 
which  can  be  picked  as  an  acceptable  and  “valid”  design.  Most  decisions  on  the 
choice  of  the  “best”  or  “suitable”  design  are  based  on  preferential  decision  making 
where  the  decision  maker  already  has  set  priorities  in  the  objective  functions.  One 
possible  approach  for  selection  of  a  “suitable”  design  could  be  done  by  probabilistic 
methods.  A  reliability  based  procedure  has  definite  advantages  over  a  deterministic 
design  procedure. 

A  general  approach  for  estimating  the  variance  of  structural  response  variables, 
given  the  mean  values  and  variances  of  system  properties  which  are  probabilistic 
in  nature,  is  outlined.  Here  the  probabilistic  properties  of  a  system  are  viewed  as 
discrete  random  variables.  The  thickness  of  a  plate  or  any  design  variable  may  vary 
due  to  mantifacturability  or  some  other  criteria,  and  hence  can  be  characterized  by 
a  mean  value  and  a  single  value  of  the  variance. 

The  statistical  parameter  of  interest  in  the  present  work  are  the  thicknesses  of 
the  plates  which  are  the  design  variables  of  the  structure.  The  structural  response 
quantities  that  we  are  interested  in  are  displacement,  stress,  and  natural  frequency. 
The  results  from  these  calculations  can  be  interpreted  in  may  different  ways.  One 
of  the  most  important  being  that  '  he  variance  calculations  generate  a  single  scalar 
quantity  for  each  of  the  different  parameters  such  as  displacement  at  a  location  of 
interest,  stress  at  a  critical  region  or  system  natural  frequency.  This  makes  it  easier 
in  judging  the  behavior  of  the  structure  to  the  variance  in  the  design  variables. 
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In  the  optimization  problem  the  constraint  function,  gi{X),  represents  each 
mode  of  malfunction  or  failure  of  the  system.  For  given  values  of  random  variables 
this  function  separates  the  space  of  random  variables  into  safe  and  failure  zones. 
That  is 

giiX)  >  0  -♦  safe  —  space  (134.a) 

gi(X)  <  0  — »  failure  —  space  (134.6) 

where  the  boimdary  between  these  two  regions,  =  0,  is  called  a  failure  surface 

or  limit  state  function.  If  the  joint  probability  density  function  of  the  random 
variables,  is  known,  then  the  reliability  measure  associated  with  the  t** 

constraint  can  be  obtained  as  follows  (Rao,  1984): 

P»i  =  [  fx{x)dx  (135) 

Jgi>o 

As  the  joint  distribution  of  the  random  variables  is  rarely  known  and  also  the  integral 
evaluation  is  extremely  diiiicult,  approximations  for  the  above  integral  have  been 
proposed  by  many  researchers.  Approximating  Eq  135,  the  probability  is  calculated 
as 

=  'Km)  (136) 


where  ^(‘j  is  the  st^mdard  normal  probability  distribution  function,  and  iji  is  the 
reliability  index  for  the  failure  mode  defined  as 


Sim 


(137) 


gii^X)  is  the  mean  value  of  the  constraint,  and  o’(^i)  is  the  standard  deviation.  These 
are  calculated  as  shown  below; 


giiX)  =  giiX) 


(138) 
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(139) 


where  n  denotes  the  total  number  of  design  variables.  The  method  for  evaluating 
7]i,  by  expanding  around  the  mean  value  of  the  random  variables,  is  known  as  a 
first-order  second  moment  method  (Nikolaidis  and  Burdisso,  1988). 

For  limit  state  functions  that  are  not  highly  nonlinear,  the  probability  of  failure 
can  be  approximated  with  good  accuracy  as 

Pfi  =  (140.a) 

or 

Pfi  =  l-  P,i  (140.6) 

The  individual  failure  probabilities  of  each  fiulure  subsystem  were  combined  to  eval¬ 
uate  the  system  probability  of  failure.  The  exact  evaluation  of  the  failure  probability 
is  a  formidable  task  that  reqiiires  knowledge  of  the  joint  probabilities  corresponding 
to  all  combinations  of  the  failure  modes. 

Approximating  the  system  failure  probability  by  various  bounds  that  can  be 
evaluated  inexpensively,  Ditlevsen’s  first  order  upper  bound  is  selected  for  calcu¬ 
lating  system  reliability  (Yang  and  Nikolaidis,  1990).  This  is  because  in  this  for¬ 
mulation  the  failure  modes  are  completely  independent,  and  the  approximation  is 
valid  for  the  range  of  values  of  Pj  usually  encountered  in  structural  optimization 
applications  (Davidson,  et.al,  1977).  System  probability  can  be  defined  as 

P,  =  UZiil  -  Pfi)  (141.tt) 
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If  the  values  of  the  Pfi  are  small,  higher  order  terms  se  neglected,  and  the  above 
equation  can  be  redefined  as  follows: 


P,  =  l 


(141.6) 


P/  =  E  Pfi 

1=1 

where  Pf  is  the  system  failure  probability,  Pfi  is  the  probability  of  the  individual 
failure  mode,  and  m  denotes  the  number  of  failure  modes. 
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Fig.  4.2.1  Multiobjective  Compound  Scaling  (MCS) 
Algorithm  Flow-Chan 
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Fig.  4.2.2  Oemonstration  problem  for  bicriterion  optimization 


✓  t 


Fig.  4.2.3  Demonstration  problem  for  bicriterion  optimization 


SECTION  5 


NUMERICAL  RESULTS  AND  DISCUSSIONS 


In  this  section  eight  example  problems  are  solved  using 


the  MCS  technique. 


5.1  Example  1:  Two  Bar  Truss 


The  two  bar  truss  shown  in  Figure  5.1.1  is  considered  for  illustrating  the  MCS 
approach  (Rao,  et.al.,  1988).  The  cross  sectional  area  of  the  member  and  the  posi¬ 
tion  of  joints  1  and  2  are  treated  as  the  design  variables.  The  truss  is  assumed  to 
be  symmetric  about  the  y  axis.  The  coordinates  of  joint  3  are  held  constant.  The 
weight  of  the  truss  and  the  displacement  of  joint  3  are  considered  as  the  objective 
functions  Fi  and  F2,  respectively.  The  stresses  induced  in  the  members  are  con¬ 
strained  to  be  less  than  the  permissible  stress,  (tq.  In  addition  ,  lower  bounds  are 
placed  on  the  design  variables.  Thus  the  problem  can  be  stated  as: 


Minimize 


Fi{X)  =  2pkx2y/l  +  XI 2 


F2{X)  = 


Ph{l+xx^y^l  +  xi*f^ 

2\'2Exi'^X2 


(142.o) 

(142.6) 
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Subject  to 


9l(X)  = 

92m  = 


f(l+xi)(H-xi»)°'° 

2\/2®i®2 

P(ll-l)(l+Xl^)°" 

2^2xiX2 


—  a-Q  <0 


-  ao  <0 


(143.a) 

(143.6) 


where  xi  =  51X2  =  E  =  Young’s  modulus,  p  —  density  of  the  material, 

the  lower  boimd  on  the  first  design  variable  is  O.ltn,  and  the  lower  bound  on  the 
second  design  variable  is  1.0  in.  The  problem  data  is  p=0.283  h  =  100  in, 
■Amin  =  1  in*,  P  =:  10,000  Ih  ,E  =  30  X  10*j^,  ctq  =  20,000  The  constraints 
are  normalized  by  dividing  the  values  with  <t^.  The  multicriteria  optinoization  was 
performed  using  the  generalized  scaling  algorithm.  The  MCS  generates  a  number  of 
optimal  Pareto  solutions.  Which  one  of  these  designs  is  to  be  taken  as  the  optimum 
is  decided  by  calculating  the  reliability  criteria  of  the  system.  After  elimination 
based  on  reliability  calctilations  of  all  the  constraints,  two  possible  designs  Table 
5.1.1,  are  obtained  which  need  some  discussion.  Case  1:  (iter  no.  7)  F\  =  86.761 
and  F2  =  0.115.  This  design  is  in  the  feasible  region  with  both  the  constraints 
being  inactive.  Case  2:  (iter  no.  8)  F\  —  81.860  and  F%  =  0.041.  This  design 
also  is  in  the  feasible  region  with  none  of  the  constraints  active.  Using  the  system 
reliability  information  the  design  at  iteration  8  was  taken  as  the  optimum.  The 
corresponding  optimum  design  variables  are  xi  =  0.858  and  X2  =  1.0972.  Optimum 
obtained  by  Rao  are  shown  in  comparison  with  the  MCS  results  in  Table  5.1.1.  For 
this  example  problem  MCS  did  not  perform  well  in  generating  some  of  the  Pareto 
optimal  solutions,  because  the  optima’s  do  not  lie  on  the  constraint  surface. 
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5.2  Example  2:  Production  Planning  Problem 

Consider  a  production  planning  problem  example  (Osyczka,A.,  1984)  in  which 
product  A  is  of  high  quality  and  product  B  is  of  lower  quality.  The  profits  are  4 
and  5  dollars  per  product  respectively.  The  best  customer  of  the  company  wishes 
to  have  as  many  as  possible  of  type  A  product.  Then  the  two  objectives,  (1)  the 
maximization  of  profit  and  (2)  the  maximum  production  of  product  A,  should  be 
considered.  Hence,  the  design  variables  are  xi  =  number  of  units  of  products  A  to 
produce  and  X2  =  number  of  units  of  product  B  to  produce.  Both  products  require 
time  in  two  departments.  Product  A  requires  1  hour  in  the  first  department  and 
1  j  hours  in  the  second  department.  Product  B  requires  1  hour  in  the  first  depart¬ 
ment  and  I  in  the  second  department.  The  available  hours  in  each  department  are 
200  monthly.  Furthermore,  there  is  a  maximum  market  potential  of  150  units  for 
product  B.  These  restrictions  are  the  constraints.  Hence,  the  optimization  problem 
can  be  written  as  follows: 

Minimize 

Fi{X)  =  1000.0  -  4x1  -  5x2  (144.a) 

F2{X)  =  1000.0  -  XI  (144.6) 

Subject  to 

^i(X)  =  xi -f-X2  <  200  (145.a) 

g2{X)  =  Ijxi  +  ^X2  <  200  (145.6) 

^3(^)  =  <  150  (145.c) 
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>  0 


(145.(i) 


In  this  example  as  the  problem  is  two-dimensional  a  graphical  representation  of 
the  problem  formulation  is  shown  in  Figure  5.2.1  (a).  In  this  space  the  inequality 
constraints  form  a  set  of  feasible  solutions.  The  Pareto  set  for  the  problem  is  shown 
in  Figure  5.2.1  (b),  from  which  it  can  be  seen  that  F**  are  the  possible  Pareto 
solutions.  The  possible  Pareto  sets  generated  by  the  MCS  are  shown  in  Table  5.2.1. 
It  can  be  seen  that  the  Pareto  solutions  generated  by  MCS  follow  the  dark  line 
shown  in  Figure  5.2.1(b)  which  are  the  Pareto  solutions  for  the  problem.  For  this 
particular  problem  the  MCS  generated  almost  the  complete  Pareto  set,  because  the 
problem  was  linear  in  nature  and  the  Pareto  solutions  lie  on  the  constraint  surfaces. 

5.3  Example  3  :  Apartment  Problem 

This  problem  is  taken  from  Eschenauer’s  et.al.,  (1990)  book.  The  problem  is 
stated  as:  maximize  floor  area  and  minimize  the  construction  cost.  The  design 
variables  are:  xj  living  room  length,  xj  hall  length,  X3  bedroom  number  three 
length,  X4  kitchen  length,  and  X5  bedroom  length. 

The  constraints  are 


4.28x1  <  30 

(146.0) 

4.28x1  >  15 

(146.5) 

3.05X4  <  12 

(146.c) 

3.05x4  >  5 

(146.d) 
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1.83x2  <  7.2 


(146. e) 


2.45x2  <  6.5  (146./) 

2.45x2  >  4.5  (146.(7) 

3.05x5  <  18  (146./i) 

3.05x5  >  10  (146.i) 

3.35x3  <  18  (146.j) 

3.35x3  >  10  (146.*) 

3.98x3  <  18  (146./) 

3.98x3  >  10  (146.m) 

*1  +  ®a  —  ®4  ~  =  0.0  (146.n) 


The  above  problem  is  a  multicriteria  linear  programming  formulation.  Table  5.3.1 
shows  the  Pareto-optimal  solutions  generated  by  MCS,  and  Figure  5.3.1  shows  the 
Pareto  solutions  generated  by  the  constraint  method.  Prom  both  of  these  we  can 
see  that  MCS  generated  a  partial  Pareto-optimal  set. 

5.4  Example  4:  Cantilever  Beam  with  an  End  Moment 

A  beam  of  10  in  length  and  1  in  (Figure  5.4.1)  wide  is  optimized  for  minimum 
weight  and  maximum  first  natural  frequency.  The  material  properties  are:  E  = 
lO'^pst,  u  =  0.3,  p  =  0.1  The  beam  is  subjected  to  an  end  moment  of  450 
lb  —  in.  The  allowable  normal  stress  is  (Tmax  —  30,000  psi.  It  is  required 

that  the  tip  deflection  be  limited  to  0.5  in.  The  lower  bound  on  the  thicknesses 
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was  0.1  in.  The  beam  is  modelled  into  20  elements  as  a  20  x  1  mesh.  In  the  initi^ll 
design  all  the  element  thicknesses  were  0.5  in.  The  initial  weight  was  0.5  lb,  and 
the  corresponding  natural  frequency  was  wi  =  999.08  For  the  optimization 
problem  the  total  number  of  design  variables  was  20  (element  thickness),  and  the 
total  niimber  of  constraints  was  82  (  42  displacements  and  40  normal  stresses). 

The  individual  optimum  obtained  for  the  two  objective  functions  are  shown 
in  the  Table  5.4.1.  Using  the  constraint  method  the  complete  Pareto  set  for  the 
problem  was  generated.  The  feztsible  designs  generated  by  the  MCS  algorithm 
and  the  Pareto  set  generated  by  the  constraint  method  are  shown  in  Figure  (5.4.2). 
From  the  figure  it  can  be  inferred  that  the  feasible  designs  generated  by  MCS  follow 
quite  closely  to  the  constraint  method.  As  the  optimization  progresses  utilizing 
the  sensitivity  tind  constraint  information  ,  system  reliability  is  calcxilated  at  each 
iteration.  Prom  Table(5.4.1)  it  can  be  seen  that  the  designs  which  need  further 
investigation  are  the  ones  from  iteration  6,  8,  12,  etc.,  since  the  reliability  of  the 
whole  structtire  is  highest  at  these  designs.  Let  us  consider  iteration  number  6. 
At  this  design  the  stress  constraint  at  the  tip  element  is  active.  The  displacement 
constraint  is  inactive.  The  optimum  obtained  for  the  two  criteria  are  W  =  0.403/6 
and  =  1123.363  For  the  design  at  iteration  8  the  same  constraint  is  active, 
but  the  displacement  at  the  tip  is  lower.  The  optimum  are  W  =  0.576  lb  and 
ufi  =  1£ 85.205  From  the  calculation  of  the  reliability  criteria  for  the  system 

we  can  select  the  design  of  iteration  6  as  it  has  a  higher  system  reliability  than  the 
other  Pareto  optimal.  Notice  for  a  9  %  increase  in  weight  a  31  %  increase  in  wi  was 
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achieved.  The  optimal  distribution  for  the  design  at  iteration  6  and  8  is  shown  in 
Table.  5.4.2. 

5.5  Example  5  :  Square  Plate  Clamped  at  the  Edges 

A  square  plate  clamped  at  the  edges  is  minimized  for  weight  and  for  maximum 
nodal  displacement  which  for  this  problem  happens  to  be  at  the  center  of  the  plate 
(node  121).  The  plate  is  subjected  to  a  uniform  lateral  load  of  110.3  kPa.  The 
plate  dimensions  are  508  x  508  mm.  The  configuration  is  shown  in  Figure  5.5.1. 
Only  a  quarter  of  the  plate  is  considered  due  to  symmetry.  The  material  pi  operties 
are  (Aluminum)  E  =  72.4  GPa,  v  =  0.3,  p  =  2.77  gm/cm}.  The  plate  was  modelled 
with  100  elements  with  a  mesh  size  of  10  x  10.  The  maximum  allowable  normal 
stress  (o'xftrji)  was  <rmax  =  82.7  MPa,  The  lower  bound  on  the  design  variable  was 
1.0  mm. 


The  initial  weight  of  the  structure  was  L.36  kg,  and  the  corresponding  displace¬ 
ment  (  at  node  121)  was  3.15  mm  with  an  uniform  design  of  7.62  mm.  The  total 
number  of  design  variables  was  100  (element  thickness),  and  the  total  number  of 
constraints  was  200  (normzd  stresses  in  the  x  and  y  direction).  No  attempt  was  made 
to  reduce  the  number  of  design  variables  by  linking.  The  MCS  technique  generated 
a  partial  Pareto  set  with  four  designs  which  need  further  consideration.  As  the  niim- 
ber  of  active  constraints  increases,  the  number  of  Pareto  solutions  decrease  which 
can  be  seen  in  this  particular  problem  wherein  the  number  of  stresses  reaching  the 
limit  is  always  ten  or  more  at  any  of  the  Pareto  optimal  designs.  These  designs  are 
shown  in  Table  5.5.1.  The  table  also  shows  the  corresponding  system  reliability  for 
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each  particular  design.  From  the  table  it  can  be  seen  that  the  design  at  iteration 
11  has  the  highest  reliability  when  compared  to  the  others.  The  partial  solution 
set  produced  by  MCS  and  the  constraint  method  are  plotted  in  Figure  5.5.2.  The 
constraint  method  shows  the  complete  Pareto  set.  At  the  optimal  design  selected 
(selected  using  reliability)  the  thickness  distribution  around  the  center  of  the  plate 
was  high.  Stresses  at  the  regions  about  the  central  axis  where  double  symmetry 
.>£  the  structure  exists  were  active.  The  thickness  distribution  at  the  optimum  is 
shown  in  Table  5.5.2.  The  optimum  design  is  plotted  in  Figure  5.5.3. 


5.6  Example  6  :  Plate  on  Four  Diagonal  Supports 

A  rectangular  steel  plate  of  10  in  x  12  in  (Figure  5.6.1)  supported  at  four  points 
on  the  diagonals  is  the  next  problem  solved.  The  plate  is  loaded  along  the  central 
and  boundary  lines  of  the  plate.  The  problem  is  simple  enough  to  analyze,  however, 
because  of  the  overhang  a  complex  saddle  deformation  pattern  forms  around  the 
support  points  (Prasad  and  Haftka,  1979).  The  material  properties  are  Yoimg’s 
Modulus,  E  =  28  X  lO^psi,  Poisson’s  ratio,  v  —  0.3  and  mass  density,  pm  =  7.71  x 
10  The  loading  configuration  is  shown  in  Table  5.6.1.  The  finite  element 

model  consists  of  49  elements  (  a  mesh  of  7  x  7  )  ,  as  only  one-quarter  of  the  plate 
is  actually  modeled  due  to  symmetry.  The  diagonal  support  point  is  assumed  to  be 
located  at  two  fifths  of  the  distance  between  the  plate’s  comer  and  its  center.  A 
total  of  6  lb  of  non- structural  mass  was  added  to  the  structure.  The  lower  bound 
on  the  design  variables  was  0.02  in.  The  limits  on  displacement  and  stress  were 
0.02  in  and  25000.0  psi,  respectively.  The  initial  design  started  from  a  feasible 
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point  and  had  a  weight  of  12.340  lb.  The  corresponding  fundamental  frequency  was 
ci;i  =  5400.662  when  all  the  design  variable  thicknesses  were  at  0.71  in. 

The  structural  weight  was  minimized  and  the  fundamental  frequency  was  max¬ 
imized  subject  to  stress  and  displacement  constraints.  The  total  number  of  design 
variables  was  49  (element  thickness),  and  the  total  number  of  constraints  was  162 
(  64  displacements  and  98  normal  stresses).  Figure  5.6.2  shows  the  Pareto  optimal 
set  generated  by  the  constraint  method,  and  the  partial  Pareto  set  generated  by 
the  MCS  algorithm.  Table  5.6.2  gives  the  individual  minima’s  and  a  few  Pareto 
optimal  with  their  respective  system  reliability.  Going  through  the  Pareto  optimal 
solution  the  design  at  iteration  no.  31  is  selected  as  the  best  compromise  design.  At 
this  design  the  displacement  at  node  8  was  active,  smd  a  few  stresses  were  active. 
None  of  the  optimal  thicknesses  were  at  the  lower  bound.  Prom  this  design  it  can  be 
seen  that  we  have  achieved  a  good  compromise  design  between  the  two  individual 
minima’s.  Figure  5.6.3  gives  the  mode  shapes  corresponding  to  the  first  and  second 
frequencies.  Table  5.6.3  gives  the  first  five  frequencies  at  the  imtial  and  optimum 
designs.  Table  5.6.4  gives  the  optimal  thickness  distribution.  Figure  5.6.4  shows 
the  optimal  design. 

5.7  Example  7:  Simply  Supported  Square  Plate 

A  square  plate  simply  supported  at  the  edges  with  dimensions  10  in  x  10  in 
was  optimized  for  miniTnuTn  weight  (Figure  5.7.1).  The  material  properties  are: 
Young’s  modulus,  E  =  lO^psi,  Poisson’s  ratio,  u  =  0.3  and  mass  density,  pm  — 

2.587x10“®^^.  A  uniformly  distributed  load  of  100  pat  was  applied.  The  plate  was 
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modelled  using  64  finite  elements  (  a  mesh  of  8  x  8  )  .  Double  symmetry  of  the  plate 
was  observed  about  the  central  axes  and  hence  only  a  quarter  model  was  considered. 
Each  element  was  controlled  by  an  independent  design  variable.  The  objective  was 
to  minimize  the  structural  weight  and  the  maximum  nodal  displacement  (node 
1)  subject  to  normal  stress  constraints.  The  limits  on  the  normal  stresses  were 
30000.0  psi.  A  uniform  plate  with  thickness  equal  to  0.5  in  was  taken  as  the  initial 
design,  which  resulted  in  an  initial  weight  of  1.25  lb  with  a  corresponding  maximum 
displacement  (node  1)  of  0.037  in  and  a  design  that  was  in  the  feasible  region.  The 
lower  bound  on  the  design  variables  was  0.1  in.  The  MCS  algorithm  generated  five 
possible  solutions.  Table  5.7.1  shows  the  five  possible  designs  *ind  their  respective 
system  reliability.  The  best  possible  design  selected  was  the  one  at  iteration  29. 
Figure  5.7.2  plots  the  MCS  algorithm  and  the  constraint  method  Pareto  solutions. 
One  will  notice  that  for  this  particular  problem  MCS  performed  extremely  well. 
The  optimal  thickness  distribution  is  shown  in  Table  5.7.2 

5.8  Example  8  :  Square  Plate  with  Internal  and  Edge  Supports 

The  final  plate  problem  was  a  square  plate  with  dimensions  12  »n  x  12  in  (Figxire 
5.8.1).  The  plate  is  supported  along  its  edges  and  also  at  internal  locations  to 
represent  a  platform  structure  with  several  internal  supports  (a  total  of  88  simply 
supported  nodes  out  of  169  nodes).  The  plate  can  be  considered  to  have  been 
made  up  of  nine  blocks.  The  material  properties  are  those  of  Aluminum  with 
Young’s  Modulus,  E  =  lO’^psi,  Poisson’s  ratio,  u  =  0.3,  and  mass  density,  pm  = 
2.587  X  10“®^^^.  The  plate  was  modelled  with  144  elements  as  a  complete  plate 
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because  of  the  inherent  properties  of  the  plate.  The  first  mode  shape  of  the  plate 
has  a  quarter  symmetry.  The  second  and  the  third  modes  have  diagonal  skewed 
synunetry  and  the  fourth  mode  has  diagonal  symmetry.  Figure  5.8.3  shows  the 
mode  shapes  corresponding  to  the  first  and  second  natural  frequencies.  A  total 
of  24  lb  of  non-structural  mass  was  added  to  the  structure  at  all  the  free  nodes. 
At  the  initial  design  all  the  variables  were  taken  as  1.0  in,  which  resulted  in  an 
initial  weight  of  38.4  lb  with  the  resulting  fundamental  frequency  u>i  =  29816.384 
The  plate  was  optimized  for  minimum  weight  and  maximum  fundamental 
frequency  subject  to  constraints  on  the  second  and  third  natural  frequency.  The 
two  constraints  were  posed  as  u)2  >  35000.0  ^  and  0/3  >  40000.0  A  minimum 
gauge  value  of  0.1  in  was  imposed  on  all  the  design  variables.  The  initial  design 
was  in  the  infeasible  region.  The  Pareto  optimal  solutions  for  the  problem  using 
the  MCS  algorithm  and  the  Constraint  method  are  plotted  in  Figure  5.8.2.  Table 
5.8.1  shows  the  possible  designs  with  their  system  reliability.  Prom  table  5.8.1  the 
best  design  selected  is  at  iteration  38.  Constraint  two  was  active  at  the  optimum. 
At  this  design  an  increase  of  1.2%  in  the  weight  brought  about  a  5  %  increase  in 
the  fundamental  frequency.  Table  5.8.2  gives  the  first  four  frequencies  at  the  initial 
and  optimum  designs.  Table  5.8.3  shows  the  optimal  thickness  distribution  values 
for  the  optimal  design.  In  general,  elements  at  the  ceni  er  of  each  block  were  thicker 
than  the  others. 
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Table  5.1.1 

Two  bar  truss  -  Pareto  optimum  set 


Iter  No. 

Reliability 

Objective  Func. 

6 

0.4996 

[68.38,  0.328] 

7 

1.0000 

[86.76,  0.115] 

8 

1.0000 

[81.86,  0.041] 

(a) 


Quantity 

Weighting 

Method 

Rao  (1987) 

Game  Theory 
Rao  (1987) 

Present 

Approach 

X«x^ 

X2 

0.7635 

1.0540 

0.7681 

1.1408 

0.8586 

1.0972 

F-F, 

Fa 

75.0595 

0.0442 

81.4137 

0.0408 

81.8609 

0.0414 

^2 

9747.6 

1306.9 

8996.0 

1180.1 

9194.0 

698.0 

Reliability 

1.0000 

1.0000 

1.0000 

(b) 
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Iter  No.  X  =  [  X  .  X  ]  F  (x)  =  t  F  (x),  F,  (x)  ]  Active  Constraint 
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Rg.  5.3.1  Pareto  optimal  solution  for  the  Apartment  problem 
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Table  5.3.1 


Apartment  Problem  -  Pareto  optimum  set 


Pareto-optimal  Solutions 


Solution  No. 

Area  [m]  ^ 

Costs  [units] 

1 

79.06 

1248.56 

2 

78.26 

1212.12 

3 

73.89 

1169.26 

4 

71.58 

1133.54 

5 

65.08 

1068.39 

6 

61.09 

1028.46 

Design  Variable  [m] 


Solution  No. 

CM 

X 

X 

3 

X  4 

^5 

1 

3.61 

2.65 

4.52 

2.11 

3.70 

2 

3.51 

2.65 

4.52 

1.74 

3.70 

3 

3.68 

1.87 

4.52 

1.68 

3.70 

4 

3.51 

1.84 

4.41 

1.64 

3.70 

5 

3.51 

1.84 

3.53 

1.64 

3.70 

6 

3.51 

1.84 

2.99 

1.64 

3.70 

^  3 


t 


S’  VI VI 
5  =  0 

{/3 


CO 

V 

2 

V 

o  ”c 

(tb  ^ 

•si  a 

•C9«-; 
«  Ft  o  o 


(«  II  II  II 
S  b]  >  Q- 


s 

Jl 


s 


•o 

s 


u 


85 


2500- 


o 

O 

o 

o 

o 

O 

c 

o 

o 

m 

o 

m 

CN 

oas/pDJ  ‘Xouanbaj^  |D|uaujDpunj 


86 


Table  5.4.1 

Pareto  optimal  solutions  tor  Cantilever  beam 


No. 

System  Reliability 

f(x)=  [f  (X^),f  (X^)] 

6 

0.0111 

[  .403,  1123.36] 

8 

0.0078 

[.576,  1885.20] 

12 

0.0078 

[.591,  1956.00] 

14 

0,0078 

[.586,  1984.00] 

26 

0.0078 

[.574,  1973.36] 

Initial 

Minimize  Wt 

Subjea  to 

Constraints 

Maximize  q) 

Subject  to  ^ 

Constraints 

Wt. 

Ob) 

0.50 

0.37 

0.78 

rad/sec 

999.08 

770.85 

2375.79 

System 

Reliability 

0.0000 

0.2086 

0.0078 
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Table  5.4.2 

Optimal  design  of  Cantilever  Beam 


Element 

Number 

Iteration  No. 

6 

Iteration  No. 

8 

1 

0.3000 

0.3000 

2 

0.3000 

0.3000 

3 

0.3000 

0.3000 

4 

0.3000 

0.3000 

5 

0.3000 

0.3000 

6 

0.3000 

0.3000 

7 

0.3027 

0.3000 

8 

0.3784 

0.4873 

9 

0.3784 

0.6092 

10 

0.4730 

0  6092 

1 1 

0.4730 

0.7614 

12 

0.4730 

0.7614 

13 

0.4730 

0.7614 

14 

0.4730 

0.7614 

15 

0.4730 

0.8184 

16 

0.4730 

0.8184 

17 

0.4730 

0.8184 

18 

0.4730 

0.8184 

19 

0.4730 

0.8184 

20 

0.4730 

0.8184 

Material  Properties: 
E  =  72.4  GPa 
=03 

p  =  277  gm/cm 
UDL=110.3kPa 


Element  Number 
X  Grid  Point  Number 


(b) 


Rg.  5.5.1  Clamped  Square  Plate  With  Uniform 
Lateral  Load 


89 


MCS _ 

Constraint  method 
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Fig.  5.5.2  Minimal  curve  in  criterion  space  (Example  5) 


Fig.  5.5.3  Optimal  Design  (Example  5) 
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Table  5.5.1 

Pareto  optimal  solutions  tor  the  Clamped  Square  Plate 


Iter  No. 

System  Reliability 

f(x)=  [f  (x),  f^Cx)] 
[Weight,  Displacement] 

11 

1.0000 

[1.2820,  1.7676] 

19 

0.9905 

[1.2426,  1.9190] 

28 

0.8518 

[1.1820,  1.9792] 

29 

0.9918 

[1.1838,  1.9717] 

Initial 

Minimize  Wt 
Subject  to 

Stress  constraints 

Minimize  maximum  displacement 
Subject  to 

Stress  constraints 

Wt. 

(Kg) 

1.36 

0.99 

11.06 

u 

(mm) 

3.1517 

2.49 

0.69 

System 

Reliability 

0.0000 

0.7628 

1.0000 
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Tabie  5.5.2 

Optimal  design  of  Clamped  square  plate 


Ele  No. 

Thickness 

1 

1.2517 

2 

1.2517 

3 

1.2517 

4 

7.5306 

5 

9.4372 

6 

8.9890 

7 

8.7843 

8 

9.1833 

9 

9.9334 

10 

9 . 1833 

11 

1.2517 

12 

1.2517 

13 

1 . 7093 

14 

1.5636 

IS 

4 . 7489 

16 

9.8852 

17 

10.6869 

18 

8.6516 

19 

7.6358 

20 

9.3592 

21 

1.2517 

22 

1.7093 

23 

1.4676 

24 

6.4603 

25 

1.5636 

26 

4.6286 

27 

4.4490 

28 

8.4229 

29 

3.3411 

30 

3.3411 

31 

7.5306 

32 

1.5636 

33 

6.4603 

34 

1 . 7722 

35 

1.4664 

36 

1.5204 

37 

1.4664 

38 

3.5196 

39 

3.3411 

40 

5.6592 

41 

9.4372 

42 

4.7489 

43 

1 . 5636 

44 

1.4664 

45 

1.4697 

46 

1.4176 

47 

1.4664 

48 

10.2041 

49 

5 . 1504 

50 

9.2838 

Ele.  No. 

Thickness 

51 

8.9890 

52 

9.8852 

53 

4.6286 

54 

1.5204 

55 

1.4176 

56 

1.0000 

57 

1 . 4664 

58 

5 . 1504 

59 

7.3205 

60 

17.4012 

61 

8 . 7843 

62 

10.6869 

63 

4 .4490 

64 

1 . 4664 

65 

1 . 4664 

66 

1.4664 

67 

4 . 4643 

68 

4 . 7489 

69 

14.9552 

70 

20.4048 

71  9.1833 

72  8.6516 

73  8.4229 

74  3.5196 

75  10.2041 

76  5.1504 

77  4 . 7489 

78  1.0000 

79  14.9552 

80  32.0786 


81 

9.9334 

82 

7.6358 

83 

3.3411 

84 

3.3411 

85 

5 . 1504 

86 

7.3205 

87 

14.9552 

88 

14.9552 

89 

1.0000 

90 

29.4741 

91 

9.1833 

92 

9.3592 

93 

3.3411 

94 

5.6592 

95 

9.2838 

96 

17.4012 

97 

20.4048 

98 

32.0786 

99 

29.4741 

100 

13.9009 
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Material  Properties 

E  =  2a.E07  psi 
1/  *  0.3 

7.71E-04  lb  s'/  in* 

Minimize  Weight 
Maximize  CJj 

Subject  to 
U  <  0.02  in 
<7  <  25000.0  psi 


Quarter  Model 


Fig.  5.6.1  Plate  Supported  Along  Four 

Diagonal  Points 
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Fig.  5.6.2  Minimal  curve  in  criterion  space  (Example  6) 


(a)  Mode  1 


(b)  Mode  2 


Fig.  5.6.3  Mode  shapes  for  Example  6 
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Table  5.6.2 


Pareto  optimal  solutions  for  the  Plate  with 
Four  Diagonal  supports 


Iter  No. 

System  Reliability 

f(x)=  [f  (X^).  f(X2)  ] 

1  9 

0.5357 

(10.34,  3651.89] 

22 

0.5337 

(10.44,  3793.82] 

25 

0.5459 

(10.50,  3897.54] 

28 

0.5394 

(10.62,  4061.11] 

31 

0.9116 

(10.26,  3725.74] 

34 

0.4959 

(10.29,  3788.21] 

36 

0.4631 

(10.40,  3923.56] 

38 

0.4615 

(10.49,  4033.87] 

Initial 

Minimize  Wl 
Subjeato 

Constraints 

Maximize  (n 

Subject  to  ^ 

Constraints 

WL 

Ob) 

“1 

(rad/sec) 

System 

Reliability 

12.35 

5400.34 

1.0000 

9.26 

1924.00 

0.3423 

14.56 

8758.00 

1.0000 

Table  5.6.3 

First  Five  Frequencies  (  Example  6) 


Initial 

Optimal 

5400.338 

3672.200 

(D 

2 

8157.842 

4420.065 

“3 

13307.113 

6790.055 

17350.319 

8775.119 

19396.273 

11241.834 
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Table  5.6.4 

Optimal  design  of  Plate  with  four  diagonal  supports 


Ele  No. 

Thickness 

Ele.  No. 

Thickness 

1 

0.5014 

25 

0 . 5034 

2 

0.5515 

26 

0 . 4932 

3 

0.5014 

27 

0.2222 

4 

0 . 5459 

28 

0 . 3949 

5 

0.5292 

29 

0.3151 

6 

0 . 5249 

30 

0.3268 

7 

0.5362 

31 

0.5456 

8 

0.7735 

32 

0.4857 

9 

0.7710 

33 

0.5081 

10 

0.7621 

34 

0.5612 

11 

0 . 6089 

35 

0.2866 

12 

0 . 6503 

36 

0.3151 

13 

0.5578 

37 

0.3279 

14 

0.5581 

38 

0.5106 

15 

0 . 5969 

39 

0.4027 

16 

0 . 7483 

40 

0.5081 

17 

0 . 8349 

41 

0.3083 

18 

0 . 6627 

42 

0.1296 

19 

0 . 4836 

43 

0.3397 

20 

0.5186 

44 

0.2075 

21 

0.5106 

45 

0.2561 

22 

0.4107 

46 

0.1050 

23 

0.3268 

47 

0.2566 

24 

0 . 8349 

48 

0.1481 

49 

0.1467 
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Material  Properties: 


E  =  1.0E07  psi 
\)  =  0.3 
p  =0.1  Ib/in 
Udl  =  100  psi 


(a) 


123456789 


81 

72 

63 

Optimization: 

54 

Minimize  Wt. 

45 

Minimize  u(l) 

36 

Subject  to 

27 

a  ^  30000.0  psi 

18 

ME?H  9  X  ^  2  Element  Number 

X  Grid  Point  Number 


(b) 


Fig.  5.7.1  Simply  Supported  Square  Plate 
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0.05  H 


Table  5.7.1 


Pareto  optimal  solutions  for  square  plate 
simply  supported  at  the  edges 


Iter  No. 

System  Reliability 

F(x)=  [F  (x).^  (X)] 
[Weight,  Displacement] 

11 

0.0247 

[0.7119,  0.1608] 

17 

0.0101 

[0.6852,  0.1674] 

23 

0.0196 

[0.6630,  0.1727] 

29 

0.0895 

[0.6401,  0.1779] 

36 

0.0211 

[0.6167,  0.1836] 

B 

Minimize  Wl 
Subject  to 

Stress  constraints 

Mininoize  Maximum  Di^lacement 
Suiqect  to 

Stress  constraints 

Wt 

m 

1.2500 

0.5915 

1.2500 

u 

(in) 

0.0372 

0.2423 

0.0372 

System 

Relialnlity 

1.0000 

0.0010 

1.0000 
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Table  5.7.2 

Optimal  design  of  simply  supported  square  plate 


Ele  No.  Thickness 

0.2870 
0.3108 
0.3182 
0.3121 
0 . 2709 
0.1933 
0.1505 
0 . 1260 
0 . 3108 
0 . 2603 
0.3121 
0.2870 
0.2353 
0.2173 
0.1751 
0.  1505 
0.3182 
0.3121 
0.2710 
0.3121 
0.1829 
0.2139 
0.2419 
0.1751 
0.3121 
0.2870 
0.3121 
0.1484 
0.2176 
0.1933 
0.2531 
0.2273 


Ele.  No. 

Thickness 

33 

0.2709 

34 

0.2353 

35 

0 . 1829 

36 

0.2176 

37 

0.1260 

38 

0.1554 

39 

0.4918 

40 

0.3061 

41 

0.1933 

42 

0.2173 

43 

0.2139 

44 

0.1933 

45 

0.1554 

46 

0.3958 

47 

0.4918 

48 

0.3627 

49 

0.1505 

50 

0.1751 

51 

0.2419 

52 

0.2531 

53 

0.4918 

54 

0.4918 

55 

0.4918 

56 

0.2531 

57 

0.1260 

58 

0.1505 

59 

0.1751 

60 

0.2273 

61 

0.3061 

62 

0.3627 

63 

0.2531 

64 

0.1333 
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Material  Properties: 

E  =  1.0E07  psi 
P  =  0.3 

p  =  3.507E-O4  Ib  s'/  in' 
Minimize  Wt 
Maximize 
Subject  to 

CJ,  ^  39000.0  rad/see 
5  40000.0  rad/sec 
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Fig  5.8.1.  SQUARE  PLATE  vaiH  INTERNAL 

AND  EDGE  SUPPORTS 
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30000  - 


Mode  1 


Mode  2 


Rg.  5.8.3  Mode  shapes  for  Example  8 
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Table  5.8.1 


Pareto  optimal  solutions  for  square  plate 
with  internal  and  edge  supports 


Iter  No. 


System  Reliability 

F(x)=  [F  (x),F2  (X)] 
[Weight,  Fundamental  Frequency] 

0.9379 

[42.88,  36461.20] 

0.9;  111 

[43.65,  36859.02] 

0.9  10 

[43.77,  37549.32] 

0.9299 

[43.82,  38906.14] 

0.0212 

[44.36,  39777.20] 

(rad/sec) 

System 

Reliability 


! 

Initial 

Minimize  Weight 
Subject  to 
Constraints  on 

Maximize  (Oi 
Subject  to 
Constraints  on 

CD  >  (0  . 

2  3 

(Og  .  ©3 

38.4 

42.32 

63.82 

29816.38 

33404.12 

43610.63 

0.0000 

0.0235 

1.0000 
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Table  5.8.2 


First  Four  natural  frequencies 


Minimze  Wt. 
Subject  to 
constraints  on 

(D  .  (0 

2  ’  3 

Maximize  «>  i 
Subject  to 
constraints  on 

CO  .  (0 

2  ’  3 

Optimum 

Design 

(0  , 

33404.116 

43610.631 

36461.208 

09 

2 

36647.414 

45184.092 

37479.610 

(0 

3 

39899.018 

48042.382 

40026.954 

0) 

4 

39975.613 

52362.123 

40539.290 
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Table  5.8.3 

Optimal  design  of  Plate  with  internal  and  edge  supports 


i.  No. 

Thickness 

Ele  No 

1 

1.5455 

40 

2 

1.3760 

50 

3 

0.0866 

51 

4 

1.5432 

52 

5 

1.3185 

53 

6 

0.5840 

54 

7 

0 . 5820 

55 

8 

1.2302 

56 

e 

0.6020 

57 

10 

1.5407 

58 

11 

1.5407 

50 

12 

0.6058 

60 

13 

1.4602 

61 

14 

1.4881 

62 

15 

0.0866 

63 

16 

1.5340 

64 

17 

1.5356 

65 

18 

1.5435 

66 

10 

1.4321 

67 

20 

1.5532 

68 

21 

1.5727 

60 

22 

1.5380 

70 

23 

1.3450 

71 

24 

1.3450 

72 

25 

1.6035 

73. 

26 

1.8253 

74 

27 

1.5105 

75 

28 

0.8078 

76 

20 

0.7034 

77 

30 

0.6520 

78 

31 

1.3526 

70 

32 

1.5543 

80 

33 

1.5407 

81 

34 

1.6788 

82 

35 

1.4304 

83 

36 

0.6036 

84 

37 

2.2266 

85 

38 

1.6850 

86 

30 

1.0868 

87 

40 

1.6083 

88 

41 

2.2006 

80 

42 

0.0300 

00 

43 

0.6553 

01 

44 

2.1413 

02 

46 

0.5522 

03 

46 

0.5600 

04 

47 

1.2521 

05 

48 

1.3505 

06 

Thickness 

Ele  No. 

Thickness 

2.1182 

07 

2.0005 

1.7228 

08 

1 . 5276 

1.0478 

90 

0.6417 

1.7475 

100 

0 . 8200 

1 . 1854 

101 

2 . 5097 

1.5057 

102 

0.8326 

1.2302 

103 

0 . 8061 

0.7162 

104 

1 . 5637 

1.5035 

105 

0 . 7680 

0.0267 

106 

0.0802 

1.7030 

107 

1 .0407 

1.2638 

108 

1.7713 

0 . 6406 

100 

0.7697 

1.7226 

110 

1 . 6664 

1.0478 

111 

1.4670 

1.0305 

112 

1.5184 

1.3807 

113 

0.0401 

0 .0080 

114 

1.2647 

0.8087 

115 

0.8038 

1.4660 

116 

1 . 1506 

0.8833 

117 

0.4036 

1.0306 

118 

1.0234 

1.0000 

110 

1.1660 

1.1324 

120 

1.2182 

0.6406 

121 

1.5805 

1.5135 

122 

1.7711 

0.8000 

123 

1.6871 

0.0013 

124 

1.3356 

1.4260 

125 

1.6350 

0.0847 

126 

1.6167 

0.0847 

127 

1.7606 

1.5178 

128 

1.6365 

0.7760 

120 

1.6328 

1.3558 

130 

1.4034 

1.5315 

131 

1.2482 

1.1100 

132 

1 .0246 

2.1140 

133 

1.7056 

1.5380 

134 

1.6306 

0.0054 

135 

1.6306 

1.8455 

136 

0.7917 

1.1002 

137 

1.1763 

1.2065 

138 

0.7624 

1.2857 

130 

0.7123 

0.7885 

140 

1.1620 

1.5548 

141 

1.4388 

1.0442 

142 

1.2011 

1.1451 

143 

1.2100 

1.1151 

144 

1.7864 
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SECTION  6 


CONCLUSIONS 


A  new  technique  for  multiobjective  optimization  was  developed.  An  effort  was 
made  to  modify  the  existing  Generalized  compound  scaling  algorithm  for  multiob¬ 
jective  optimization  techniques.  This  technique  was  demonstrated  on  lai'ge  scale 
problems  with  more  than  one  hundred  design  variables  and  above  three  hundred 
constraints.  Objective  functions  which  were  linear  as  well  as  non-linear  were  used 
for  demonstration  of  the  multiobjective  algorithm.  Complete  Pareto  solution  sets 
were  generated  for  the  problems  using  the  constraint  method  for  comparison  with 
the  present  approach.  In  most  of  the  problems  the  MCS  Pareto  set  followed  quite 
closely  to  the  constraint  method.  The  technique  proved  to  be  robust  as  well  as 
computationally  efficient.  The  best  compromise  design  was  selected  based  on  the 
system  reliability  of  the  structure.  This  proved  to  be  an  efficient  way  for  selecting 
the  best  possible  compromise  solution.  As  in  most  structural  optimization  problems 
the  optimum  lies  on  the  constraint  boimdary  of  the  active  constraint  ol  at  the  inter¬ 
section  of  the  active  constraints.  This  assumption  was  the  basis  in  developing  the 
algorithm.  Pareto  optimal  solutions  also  lie  on  the  constraint  surface,  and  since  the 
objectives  are  also  treated  as  constraints,  generating  Pareto  solutions  was  straight 
forward. 
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Some  of  the  drawbacks  of  this  method  are  that  Pareto  solutions  are  found 
only  when  the  solutions  lie  on  the  constraint  boundary.  If  this  is  not  the  case, 
then  arbitrary  solutions  are  generated  which  may  or  may  not  be  Pareto  optimal. 
For  some  of  the  problems  where  one  of  the  objective  functions  is  highly  non-linear 
(frequency),  the  MCS  generated  solutions  crossed  over  the  constraint  method  curves 
hence  generating  solutions  which  were  not  true  Pareto  optimum.  This  was  seen 
for  Example  8  where  the  number  of  design  variables  was  144.  This  was  not  seen 
for  cases  when  weight  and  displacements  were  used  as  the  conflicting  objective 
functions.  On  the  whole  when  comparing  MCS  to  the  methods  described  in  Section 
1  (eg.,  constraint,  weighting,  NISE,  global  criterion  methods)  ,  the  algorithm  proved 
to  be  reliable  and  robiist  requiring  less  computer  resources  to  generate  a  Pareto 
set.  \'alidation  of  this  technique  on  plate  structures  with  displacement,  stress, 
and  frequency  ftmctions  (constraints  or  objectives)  where  the  problems  are  highly 
nonlinear  demonstrated  the  algorithms  capabilities. 


REFERENCES 


Balachandran,  M.,  and  Gero,  J.S.,  1987:  Use  of  knowledge  in  selection  and 
control  of  optimization  algorithm.  Enginee  ing  Optimization,  Vol.  12,  No.  2, 
163-173. 

Bartel,  D.L.,  and  Marks,  R.W.,  1974:  The  optimum  design  of  mechanical  sys¬ 
tems  with  competing  design  objectives.  Journal  of  Engineering  for  Industry, 
Vol.  96,  171-178. 

Belytschko,  T.,  Tsay,  C.S.,  and  Liu,  W.K.,  1981:  A  stabilization  matrix  for  the 
bilinear  Mindlin  plate  element.  Computational  Methods  in  Applied  Mechanics 
and  Engineering,  Vol.  29,  No. 3,  313-327. 

Belytschko,  T.,  and  Tsay,  C.S,,  1983:  A  stabilization  procedure  for  the  quadri¬ 
lateral  plate  element  with  one-point  quadrature.  International  Journal  for  Nu¬ 
merical  Methods  in  Engineering,  Vol.  19,  No. 3,  405-419. 

Brockman,  R.A.,  1987:  Dynamics  of  the  bilinear  Mindlin  plate  element.  Interna¬ 
tional  Journal  for  Numerical  Methods  in  Engineering,  Vol.24,  No.l2,  2343-2356. 

Brockman,  R.A.,  Lung,  F.Y.,  and  Braisted  W.R.,  1989:  Probabilistic  finite 
element  analysis  of  dynamic  structural  response.  WPAFB,  Technical  Report, 


114 


AFWAL-TR-88-2149. 


Cohon,  J.L.,  1978:  Multiobjective  programming  and  planning.  Academic  Press, 
New  York. 

Cook,  R.D.,  Malkus,  D.S.,  and  Plesha,  M.E.,  1989:  Concepts  and  applications 
of  finite  element  analysis  (Third  Edition).  John  Wiley  and  Sons,  190-192. 

Davidson,  J.W,  Felton,  L.P,  and  Hut,  G.C.,  1977:  Optimum  design  of  struc¬ 
tures  with  random  parameters.  Computers  and  Structures,  Vol.  7,  4^1-486. 

Diaz,  A.,  1987:  Interactive  solution  to  multiobjective  optimization  problem. 
IntemationalJoumcd  for  Numerical  Methods  in  Engineering,  Vol.  24,  1865-1877. 

Dhingra,  A.K.,  Rao,  S.S.,  and  Miura,  H.,  1990:  Multiobjective  decision  making 
in  a  fuzzy  environment  with  applications  to  helicopter  design.  AIAA  Journal, 
Vol  28,  709-710. 

Duckstein,  L.,  1984:  M  ultiob jective  optimization  in  structural  design:  The 
model  choice  problem.  In  E.  A  trek  et.al.  (eds.).  New  Directions  in  Structural 
Design,  WHey,  New  York. 

El-Sayed,  M.E.M,  Ridgely,  B.J.,  and  Sandgren,  E.,  1989:  Nonlinear  structural 
optimization  using  goal  prograr  iming.  Computers  oTid  Structures,  Vol.  92,  No.l, 
69-79. 


115 


Eschenauer,  H.,  Koski,  J.  and  Osyczka,  A.,  1990:  Multicriteria  design  optimiza¬ 
tion,  Procedures  and  Applications.  Springer,  Verlag,  Berlin,  Hedelberg. 

Grandhi,  R.V.,  Venugopal,  N.S.,  and  Venkayya,  V.B.,  1991:  Optimality  criteria 
method  for  minimum  weight  design  of  plate  structures.  S2nd  AIAA  /ASMS 
/ASCE  /AES  /ASC  Structures,  Structural  Dynamics,  and  Material  Conference, 
Baltimore,  MD,  Part  1,  AIAA-91-092S,  21S-226. 

Hajela,  P.  and  Shih,  C.J.,  1990:  Multiobjective  optimum  design  in  mixed  integer 
and  discrete  design  variable  problem.  AIAA  Journal,  Vol.  28,  670-675. 

Hughes,  T.J.R,  Cohen,  M.  and  Haroun,  M.,  1978:  Reduced  and  selective  in¬ 
tegration  techniques  in  finite  element  analysis  of  plates.  Nuclear  Engineering 
Design,  Vol.  20S-222. 

Hwang,  C.L.  and  Masud,  A.S.M.,  1979:  Mtiltiple  objective  decision  making  - 
Methods  and  Applications.  Springer,  Verlag,  Berlin. 

Kosld,  J.  and  Silvennoinen,  R.,  1982:  Pareto  optima  of  isostatic  trusses.  Com¬ 
putational  Methods  in  Applied  Mechanical  Engineering,  Vol.  SI,  No.  S,  265-279. 

Kosld,  J.,  1984:  Bicriterion  optimum  design  method  for  elastic  trusses.  Acta 
Polytechnica  Scandinavica,  Mechanical  Engineering  Series  No.  86,  Helsinki. 

Kosld,  J.,  1987:  Norm  methods  and  partial  weighting  in  multicriterion  optimiza¬ 
tion  of  structures.  International  Journal  for  Numerical  Methods  in  Engineering, 


116 


Voi  24,  No.  21,  1101-1121. 


Koski,  J.,  1981:  Multicriteria  optimization  in  structural  design.  International 
Symposium  on  Optimum  Structural  Design:  The  11th  ONR  Naval  Structural  Me¬ 
chanics  Symposium,  Tucson,  Arizona. 

Koski,  J.,  1984:  Multicriterion  optimization  in  structural  design.  In  E.  Atrek 
et.al.  (eds.),  Neva  Directions  in  Structural  Design,  Wiley,  Neva  York,  4SS-50S 

Kuhn,  H.W.  and  Tucker  A.W.,  1951:  Nonlinear  programming.  In  Proceedings 
of  the  Second  Berkeley  S  imposium  on  Mathematical  Statistics  and  Probability. 
(Edited  by  J.  Neyman),  4S1-491. 

Mindlin,  R.D.,  1951:  Influence  of  rotary  inertia  and  shear  on  the  bending  of 
elastic  plates.  Journal  of  Applied  Mechanics,  Vol.  18,  1031-1036. 

Nikolaidis,  E.,  and  Burdisso,  R.,  1988:  Reliability  based  optimization:  A  safety 
index  approach.  Computers  end  Structures,  Vol.  28,  No.  6,  781-788. 

Osyezka,  A.,  1984:  Multicriterion  optimization  on  engineering,  Ellis  Horwood 
Series  in  Engineering  Sciences. 

Prasad,  B.,  and  Haftka,  R.T.,  1979:  Optimum  structural  design  with  plate 
bending  flnite  elements.  ASCE  J.  Structures  Division,  Vol.  105,  STll,  2367- 
2382. 


117 


Radford,  A.D.,  Gero,  J.S.,  Rosenman,  M.A.,  and  Balachandran,  M.,  1985: 
Pareto  Optimization  as  a  Computer-Aided  design  tool.  Gero,  J.S.  (ed):  Op¬ 
timization  in  Computer-Aided  Design.  Amsterdam:  North-Holland,  .{7-69. 


Rao,  S.S.,  1984:  Multiobjective  optimization  in  structural  design  with  uncertain 
parameters  and  stochastic  processes.  AIAA  Journal,  Vol.  22,  No.l,  1670-1678. 

Rao,  S.S.,  1984:  Optimization,  Theory  and  Applications  (Second  Edition).  John 
WUey  and  Sons,  610-612. 

Rao.  S.S.,  1987:  Game  theory  approach  for  multiobjective  structural  optimiza¬ 
tion.  Computers  and  Structures,  Vol  25,  119-127. 

Rao,  S.S.,  Venkayya,  V.B.,  and  Khot,  N.S.,  1988:  Optimization  of  actively 
controlled  structures  using  goal  programming  techniques.  International  Journal 
for  Numerical  Methods  in  Engineering,  Vol  26,  No.  1,  183-197. 

Rao,  S.S.,  Venkayya,  V.B.,  and  Khot,  N.S.,  1988:  Game  Theory  approach  for 
the  integrated  design  of  structures  and  controls.  AIAA  Journal,  Vol.  26,  No. 
4,  4<iS-469. 

Rao,  S.S.,  Design  of  vibration  isolation  systems  using  multiobjective  optimiza¬ 
tion  techniques.  ASME  paper,  No.  84-DET-6O. 

Sandgren,  E.,  1989:  Structural  design  optimization  for  latitude  by  nonlinear 
goal  programming.  Computers  and  Structures,  Vol  33,  No.  6,  1395-1402. 


118 


Saravanos,  D.A.,  and  Chaims,  C.C.,  1990:  Multiobjective  shape  and  mate¬ 
rial  optimization  of  composite  structures  including  damping.  Presented  at  Slat 
Structures,  Structural  Dynamics,  and  Materials  Conference,  Long  Beach,  CA,  2-4 
Apr,  S71-S77. 

Schmit,  L.A.,  and  Farshi,  B.,  1974:  Some  approximation  concepts  for  structural 
synthesis.  AIAA  Journal,  Vol.  12,  No.  5,  692-699. 

Stadler,  W.,  1984:  Multicriteria  Optimization  in  Mechanics  (  A  Survey).  Ap¬ 
plied  Mechanics  Reviews,  Vol.  37,  No.  3,  277-286. 

Stadler,  W..  1987:  Update  to  multicriteria  Optimization  in  Mechanics  (  A  Sur¬ 
vey).  Applied  Mechanics  Update,  ASME,  4i'^~4^0. 

Starnes,  J.H,  and  Haftka,  R.T.,  1979:  Preliminary  design  of  composite  wings 
for  buckling,  stress  and  displacement  constraints.  Journal  of  Aircraft,  Vol.  16, 
564-570. 

Tseng,  C.H.,  and  Lu,  T.W.,  1990:  Minimax  mtJtiobjective  optimization  in 
structural  design.  International  Journal  for  Numerical  Methods  in  Engineering, 
Vol.  30,  No.  6,  1213-1228. 

Venkayya,  V.B.,  1971:  Design  of  optimum  structures.  Computers  and  Struc¬ 
tures,  Vol.  1,  265-309. 


119 


Venkayya,  V.B.,  1986;  A  further  Generalization  of  Optimality  Criteria  Ap¬ 
proaches  in  Structural  Optimization.  Proceedings  of  the  International  Confer¬ 
ence  on  Computational  Engineering  Science,  April  IO-I4,  Atlanta,  Ga.,  46i.l-15. 

Venkayya,  V.B.,  and  Tischler,  V.A.,  1989:  A  compound  scaling  algorithm  for 
mathematical  optimization.  WPAFB  Technical  Report,  WDRC-TR-89-3040. 


Venkayya,  V.B.,  1989:  Optimality  Criteria:  A  basis  for  multidisciplinary  design 
optimization.  Journal  of  Computational  Mechanics,  Vol.  5,  1-21. 

Venkayya,  V.B.,  Kolonay,  R.M,  Tischler,  V.A,  and  Canfield,  R.A.,  1990:  A 
generalized  optimality  criteria  for  mathematical  optimization.  A  paper  presented 
at  AIAA  Dynamics  Specialist  Conference,  Long  Beach,  CA,  April. 

Venkayya,  V.B.,  1991:  Design  optimization  in  dynamic  enviromnent.  Intenui- 
tional  Conference  on  Recent  Advances  in  Structural  Dynamics,  Southhampton, 
England,  July. 

Whitney,  J.M.,  1973:  Shear  correction  factors  for  orthotropic  laminates  under 
static  load.  Journal  of  Applied  Mechanics,  Vol.  40,  S02-  SO4. 

Yang,  J.S.,  and  Nikolaidis,  E.,  1990:  Design  of  aircraft  wings  subjected  to  gust 
loads:  A  safety  index  based  approach.  AIAA,  Vol  29,  No.5,  804-812. 


120 


11. S.  Government  Printing  Office  7Sn-113 


